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ABSTRACT 

The fundamental plane of black hole activity is a relation between X-ray luminosity, 
radio luminosity, and black hole mass for hard state Galactic black holes and their 
supermassive analogs. The fundamental plane suggests that, at low-accretion rates, the 
physical processes regulating the conversion of an accretion flow into radiative energy 
could be universal across the entire black hole mass scale. However, there is still a 
need to further refine the fundamental plane in order to better discern the radiative 
processes and their geometry very close to the black hole, in particular the source of 
hard X-rays. Further refinement is necessary because error bars on the best-fit slopes 
of the fundamental plane are generally large, and also the inferred coefficients can be 
sensitive to the adopted sample of black holes. In this work, we regress the fundamental 
plane with a Bayesian technique. Our approach shows that sub-Eddington black holes 
emit X-ray emission that is predominantly optically thin synchrotron radiation from 
the jet, provided that their radio spectra are flat or inverted. X-ray emission dominated 
by very radiatively inefficient accretion flows are excluded at the >3(T level. We also 
show that it is difficult to place FR I galaxies onto the fundamental plane because their 
X-ray jet emission is highly affected by synchrotron cooling. On the other hand, BL Lac 
objects (i.e., relativistically beamed sub-Eddington AGN) fit onto the fundamental 
plane. Including a uniform subset of high-energy peaked BL Lac objects from the 
Sloan Digital Sky Survey, we find sub-Eddington black holes with flat/inverted radio 
spectra follow log La^ = (1.45 ± 0.04) logL^ - (0.88 ± 0.06) log A/bh - 6.07± 1.10, with 
fint = 0.07 ± 0.05 dex. Finally, we discuss how the effects of synchrotron cooling of jet 
emission from the highest black hole masses can bias fundamental plane regressions, 
perhaps leading to incorrect inferences on X-ray radiation mechanisms. 

Key words: accretion, accretion discs - black hole physics - radiation 
mechanismsmon-thermal - methodsistatistical - galaxies:jets - X-rays:binaries 



1 INTRODUCTION 

Accreting supermassive black holes (SMBHs), i.e.. Active 
Galactic Nuclei (AGN), are thought to have important cos- 
mological ramifications. The immense energy output by an 
AGN can feedback with the environment near a SMBH 
and potentially influence the formation of galaxies (e.g., 
Idi Matteo et all l2005l : ICroton et al] l2006l : iFontanot et all 
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l2006l : lHardcastle et al.ll2007l '). Thus, accreting black holes are 
more than interesting laboratories to explore physics around 
extreme gravitational potentials. They are also critical com- 
ponents for any self consistent view of galaxy formation. 
However, the underlying physics of black hole accretion, 
and especially the mechanisms connecting the accretion flow 
with large scale outflows (i.e., the so-called disk/jet connec- 
tion), are still outstanding problems in astrophysics. 

An increasing number of studies over the past few years 
suggest that accretion physics scales globally across the en- 
tire black hole mass scale, from ~10 Mq Galactic black holes 
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(GBHs) to 10^-1 " SMBHs (e g.. iDone fc Gierliriskil 
20051: 1 Jested l2005l; iNipoti et al.ll2005l: iKording et al.ll2006bl: 



McHardv eral] l2006l: iKording et al. I l2007l : iMarkoff et all 



2008l : lGliozzi et al.ll201Gl : lKellv et al.ll201lf ). The discovery of 



the fundamental plane of black hole activity (hereafter the 
FP) - a relation between X-ray luminosity {Lx), radio lumi- 
nosity {La), and black hole mass (Mbh) for low-accretion 
rate black holes - is one of the most prominent pieces of 
evidence supporting the unification of GBHs an d SM BHs 
jMerloni et aLlbOOa iFalcke et aLll2004l . hereafter lM03l and 
IF04| . respectively V That is, once normalized by black hole 
mass, the process that regulates the fraction of energy ac- 
creted onto a black hole that is ultimately converted into ra- 
diation could be universal for weakly accreting black holes. 
The FP is a natural consequence if black hole accretion and 
relativistic jet physics are scale invariant. In other words, 
any spatial dependence of physical properties like magnetic 
field strength, jet power, etc., can be scaled to a charac- 
teristic size (e. g., the gravitational radius Tq — GM/c?) for 
all black holes ([Falcke fc BiermannlllQQSl : iHeinz fc SunvaevI 



l2003l : lMarkoflF et al.ll20oi r We define the FP in this work as 



log Lx = in log Lr + log Mbh + B, 



(1) 



and we refer to and as the FP coefficients. 

The FP implies that, under certain conditions, it is rea- 
sonable to apply our general knowledge of GBHs toward 
AGN, and vice versa. Opposing ends of the mass scale offer 
complementary insight, and the FP therefore allows a more 
holis tic view of b lack hole accretion and jet physics (see, 
e.g.. lMarko^l20ld ). For example, GBHs present the oppor- 
tunity to follow dramatic changes in accretion processes in 
real time: GBH outbursts and subsequent dimming back to 
quiescence typically last only months to years. During an 
outburst, we often witness the launching and quenching of 
a radio jet. GBHs also undergo state transitions marked by 
different contributions of hard non-thermal X-ray emission 
(from a jet/corona) to the total observed X-ray flux, where 
softer X-rays are emitted by the accret i on disk (see, e.g.. 
Fender fc Bellonij l2004l: lMcClintockll2004l : Ih Oman fc BeUonil 



20051 : iFender et al.ll2009^ 



The "human timescale" is a major advantage that stel- 
lar mass black holes hold over SMBHs, since any equivalent 
AGN episode would last 10^-10* times longer. On the other 
hand, we only kno w of around two to three dozen GBHs 
in our Galaxy (see lDunn et al.ll201Gl ). precluding large sta- 
tistical studies, and we cannot yet predict when GBHs will 
undergo outbursts. AGN are much more common, with over 
a million so far catalogued (|Richards et al.ll2009l ), offering 
the opportunity for statistical studies to obtain more ro- 
bust constraints on general trends. AGN central black hole 
masses also cover almost four decades in mass, thus provid- 
ing a larger dynamic range than GBHs when searching for 
empirical constraints on black hole mass scalings. 

The FP can be exploited to distinguish the radiative 
processes responsible for emission originating very close to 
an accreting black hole. The nature of this emission espe- 
cially in the X-ray waveband is sti ll a controversial sub- 
ject (see, e.g., §3.2 of iNaravanlboOSi ). Observed X-rays are 
likely a superposition of several components, potentially in- 
cluding synchrotron radiation from a jet, synchrotron self 
Compton (SSC), emission from the accretion flow, and in- 
verse Compton scattering of lower energy photons off a 



corona. All of these components likely contribute to the ob- 
served X-ray spectra at some level. However, we may ex- 
pect a speciflc co mponent(s) to dominate under certain con- 
ditions (see, e.g.. iMarkoff et al.1 [ioOll: iMarkoff et al.1 l2003l : 



iMarkoff etallbOOa " also see lGilfanov|[2oTor for a recent re- 
view on X-ray emission from X-ray binary systems). IMOSI 
show that different FP coefficients are predicted in the case 
that X-rays are predominantly optically thin synchrotron 
jet emission, compared to the scenario where X-rays origi- 
nate by inverse Compton scattering. The slope of the plane 
can further diagnose the radiative efficiency of the accre- 
tion flow. Several studies over the past few years (includ- 
ing the original discovery papers) have indeed used the FP 
to i nvestigate the origin of X-rays , but with varied suc- 
cess (IKording et al.ll2006al. hereafter IK06|: IWang et al. | [2006l : 



iLi et al.ll2008l : iGiihekin et aLllioOQal : lYuan et al.ll2009l) 



A major challenge to realizing the plane's full poten- 
tial to diagnose X-ray emission is that large (non-isotropic) 
observational uncertainties are associated with Lx, Lr, and 
especially Mbh . Thus, inferring reliable coefficients becomes 
a challenging multivariate regression problem requiring ad- 
vanced statistical methods. iKOg investigate this issue by 
exploring different regression techniques while incorporat- 
ing detailed error budgets for multiple samples of accreting 
black holes. They show that improper accounting of obser- 
vational error bars will bias the best-fit values of and 

IkogI also show that the inferred coefficients and the in- 
trinsic scatter around the FP depend on the sample used to 
fit the plane. A global radio/X-ray correlation is observed 
for quiescent G BHs in the hard (i.e., low-accre tion rate, non- 
thermal) state (jCorbel et a"l]|2003l : iGallo et al 2003) . When 
more luminous (soft-state, thermally dominated) GBHs 
are included with higher accretion-rates, th e radio jet is 
quenched and the cor relation breaks down (jFender et al.l 
ll999l : lGallo et al.ll2003l l. The FP is essentially an extension 
of the hard-state GBH radio/X-ray correlation to include a 
mass term. iKoel show that the intrinsic scatter about the 
FP is smallest when only considering GBHs and AGN with 
the lowest Eddington normalized luminosities. When adding 
higher- luminosity sources, the FP coefficients change and 
the intrinsic scatter increases. The apparent sensitivity of 
the FP on accretion rate supports the notio n that AGN 
display accretion "stat es" similar to GBHs fMcio? l200l|; 
Maccaronc ct a l. 2003: Ferider fc BcUoni 2004; Jester 20ol; 
Markowitz fc UttlevH2005l: IKording et alj|2006bl. I2007l l 



iGiiltekin et al.l (|2009al 'l further explored the sensitivity 
of the FP regression on the adopted black hole sample. Pre- 
vious FP studies used highly heterogeneous samples, partic- 
ularly for the AGN, assembled from different surveys using 
different telescopes. IGiiltekin et al.l l|2009al ') overcame such 
problems by using a (quasi)-volume-limited sample of AGN, 
with all X-ray data obtained by Chandra and analyzed uni- 
formly. They also only cons ider AGN with d ynamic al black 
hole mass estimates from IGiiltekin et al.l ( 2009lj ). which 
are more reliable t han secondary scal i ngs like reverberatio n 
mapping (se e. e.g.JOnken et"al]|2004l : I Peterson et aLlbOOj ). 
The work of IGiiltekin et al.l (|2009a[ ) is a powerful proof of 
concept: they show a carefully-crafted uniform sample has 
the potential to sharpen the plane enough to eventually yield 
physically meaningful constraints on black hole accretion 
physics and unification. However, even though their AGN 
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sample covers a large range of black hole masses and accre- 
tion rates, it only includes 18 AGN. This is much too small to 
distinguish radiation mechanisms. Interestingly, even with 
a small sample, they already confirm that the inclusion of 
high-accretion rates sources (especially Seyfert galaxies) in- 
creases the intrinsic scatter about the FP. 

There is an overwhelming need to further refine the FP. 
We must have more confidence in the derived best-fit co- 
efficients in order to place more meaningful constraints on 
the physical processes connecting accretion inflows and out- 
flows. More reliable coefficients are also required if the FP 
is to ever become a viable method to accurately estimate 
black hole masses (assuming one has radio and X-ray lu- 
minosity measurements). Here, we focus on sub-Eddington 
AGN (rh < 0.01 — 0.02) that we expect to be in an analo- 
gous spectral state to hard state GBHs. Since these objects 
display a narrower intrinsic scatter about the plane, such a 
sample will allow a more rigorous exploration on how statis- 
tics and sample selection might bias the inferred slopes of 
the FP. 

A brief summary of the theoretical derivation of the FP 
is outlined in ^ Then, for the first time, we fit the FP with 
a Bayesian multivariate regression technique. We find that 
our adopted statistical approach can more tightly constrain 
X-ray radiative processes than previous studies (33 . Then, 
in the context that X-rays are predominantly optically thin 
synchrotron jet emission, we explore the efl'ects of sample 
selection by attempting to place a large uniform sample of 
BL Lac objects on the FP {^^. We consider BL Lac objects 
because, due to beaming effects, we know that their spectral 
energy distributions (SEDs) are dominated by jet emission 
(even in the X-ray). Also, they host the most massive black 
holes (Mbh ~ 10* — 10^ Mq), so they provide the largest 
possible dynamic range when compared to GBHs. Finally, 
we compare to previous FP studies in ^ and we summarize 
our results in ^ Throughout we define spectral indices as 
fi, ~ we use M to express accretion rates in physi- 

cal units, and normalized (i.e., unitless) accretion rates are 
written as rh = M/MEdd- We adopt a standard cosmology: 
Ho = 71 km s"^ Mpc"^; SI™ = 0.27; and D,a = 0.73. 



2 DERIVING THE PLANE: THEORETICAL 
SCALINGS OF MULTIWAVELENGTH 
EMISSION WITH BLACK HOLE MASS 

The FP was independently discovered by IMOSI and IF04 
who derived the expected coefficients ^u, ^m, and B (Equa- 
tion [l]) under the assumption that the radio emission is 
synchrotron from a scale-invariant conical jetQ We refer 
the reader to the discovery papers for details. Briefly, the 
amount of synchrotron radiation emitted from a scale invari- 
ant jet depends on the black hole mass and accretion rate 



^ The mathematical representation for a scale invariant jet states 
that any quantity / required to describe the jet structure (e.g., 
magnetic field strength) can be independently separated into a 
normalization set only by boundary conditions, (f>f, and a struc- 
ture function describing how the jet changes with scaled radius, 
ipfir /rp): f{r,M,rn) = (f>f{M,m)'ipf{r/rg) jHeinz fc SunvaevI 
l2003ll . 



iFalcke fc BiermannI [l995l : iHeinz fc SunvaevI 120031 '). Substi- 
tuting the appropriate expression for accretion rate as a 
function of X-ray luminosity (which can also depend on 
black hole mass) into the expression for radio luminosity 
leads to the FP. The exact expression to use for X-ray lumi- 
nosity depends on if the X-rays are modeled predominantly 
as optically thin jet synchrotron or inverse Compton emis- 
sion, and different coefficients are predicted in each case. 

We stress that jet synchrotron and inverse Compton 
emission are not mutually exclusive. Even if a jet's opti- 
cally thin synchrotron radiation extends into the X-ray, then 
the underlying accretion flow can also still produce X-rays. 
Throughout this text we w ill make comparisons to the theo- 
retical predictions in lM03i and IF04| . and we will make state- 
ments on whether the X-rays are emitted primarily by the 
'jet' or by the 'corona.' The terms 'jet' and 'corona' are 
meant only to designate the mechanism emitting most of 
the X-rays (i.e., optically thin jet synchrotron from a mag- 
netized coUimated outflow vs. inverse Compton off a corona, 
respectively) , and should not be misinterpreted to mean only 
a jet or only a corona is present. In fact, the presence of mul- 
tiple components contributing to the observed X-rays likely 
contributes to the observed intrinsic scatter about the FP. 

The results of the two FP discovery papers are consis- 
tent with each other, and the two papers generally make 
similar assumptions. For example, both papers assume that 
effects due to different black hole spins from source to source 
are of second order and will be manifested only by introduc- 
ing excess scatter in the FP. Both papers also assume that 
radiative losses from synchrotron cooling are not important, 
because radiative cooling affects the emitting particles' dy- 
namics. The latter is an important assumption, as the pre- 
dicted FP coefficients for the 'jet' model are only applica- 
ble if X-rays are strictly optically thin (i.e., uncooled) syn- 
chrotron. A derivation of the FP that includ e s syn chrotron 
radiative cooling losses is presented in iHeinj (|2004l '). 

Despite the above similaritie s, th e two groups took dif- 
ferent approaches. For example, IF04| derive the FP specif- 
ically for the case of sub-Eddington black holes assuming 
X-ra ys are optically thin synchrotron emission from jets. 
IM03| consider all accretion rates, and they also derive more 
general expressions applicable not just to optically thin syn- 
chrotron X-rays, but also where X-rays are inverse Compton. 
An important difference between the two discovery papers is 
they use different notation: IM03| deflne log L_r as the depen- 
dent variabl e and flt log Lr = S^rx log Lx + (,rm log M -I- c, 
compared to|FOJ who cast log Lx as the dependent variable. 
We adopt IF04| 's no tation in this worl|3 (see Equation [l|, 
which is related to IM03| 's notation by ^b. = ^/^rx and 
£m_= —£,rm/(,rx- In Sections 12.11 and 12.21 we transcribe 
lM03l 's expressions for the FP coefficients (their Equations 
11 and 13) to our notation for both the optically thin syn- 
chrotron jet and inverse Compton models for clarity, as we 
will eventually compare our results to their theoretical pre- 
dictions. Note, the theoretical formulae for the FP coeffi- 
cients do not account for intrinsic scatter about the plane, 
so it is possible to convert between the two notations. How- 



We us e the notation of |F04| here primarily because it was also 
used by IKOGI. W e will eventually compare coefficients to those 
derived by iKOa . so this is the more natural notation for us. 
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ever, in reality there is intrinsic scatter, so one cannot regress 
a sample of black holes in one notation, and then convert 
the best-fit regression coefficients to the other notation using 
the above formulae. Instead, the regression must be repeated 
adopting a different dependent variable. 

Finally, we note that the FP has been criticized as a 
spurious correlation driven by a si milar dependence of Lr 
and Lx on distance (|Bregmanll2005l ) . We show in ^that the 
FP coefficients for sub-Eddi ngton bl ack h oles match those 
predicted by the jet model in lMOSi and lF04l . On a qualitative 
level, it would be a remarkable coincidence for a common 
dependence on distance to exactly p roduce the predicted 
non-linear correlation. Quantitatively, I Merloni et al.l (|2006l ) 
show in great detail through multiple statistical tests that 
the FP is not simply an artifact of plotting dista nce vs. dis- 
tance; other auth ors reach the same conclusion (|M03| : IK06| : 
IWang et al.ll2006D . 



2.1 Expected FP Coefficients for Optically Thin 
Synchrotron X-rays 

In the case that X-rays are dominated by optically thin syn- 
chrotron from a jet: 



{p + 4){p + 5) 



jet 



2(2p -f 13 aRp + Qur) 



'jet 



= 3 



(pH-5)(2p + 13 + 2Qfl) 
2(2p-Kl3 + aflp + 6ai?) 



(2) 



where p is the power law index of the accelerated relativistic 
electrons (i.e., drie/d-^ ~ 7~^)i with values p ~ 2 — 3 typ- 
ical, and we can substitute p = 2ax + 1 for optically thin 
synchrotron, where ax is the observed X-ray spectral index. 
Note, our lack of knowledge about the physical conditions 
in the jet is absorbed into two observable quantities, the ra- 
dio and X-ray spectral indices aR and ax- Since the radio 
and X-rays are emitted by the same underlying source, we 
expect less intrinsic scatter about the FP if X-rays are op- 
tically thin jet synchrotron compared to inverse Compton. 



2.2 Expected FP Coefficients for Inverse 
Compton X-rays 

For the case that the X-rays originate as inverse Compton 
from a corona: 



91n<I>s 
9 In 771 



g(p-f 4) 



2p + 13 + aRp + 6a J?. 



91nm 



n 



91n<l>s 
9 In A/ 



c'ln^s 
9 In 771 



g(2p-H3-H2Qfl) 
2p + 13 -I- aRp -\- QaR 



(3) 



where {M, rh) describes the boundary conditions of the 
magnetic field at the base of the jet, and q parameterizes 
the radiative efficiency of the accretion flow {Lx ~ Mm''). 
Values of g, dln^s /d\n M , and 9 In $fl/9 In tti f or diff erent 
radiation mechanisms can be found in Table 3 of lMOSl . 



2.2.1 Radiatively Inefficient Accretion Flows 

Sub-Eddington black holes are most likely fed by radiatively 
inefficient accretion flows (RIAFs), so we will primarily con- 
sider values of g > 2 when comparing FP regressions to the- 
oretical predictions (radiatively efficient accretion disks have 
q — 1). Any such mechanica lly cooled accret i on flo w should 
ha ve oc M~^rn (Heinz fc Sunvaevll20oj : iMOBj ). Follow- 
ing |mO^, the specific RIAF model we will con sider here is an 
adve c tion dominated acc retion flow (ADAF; iNaravan fc Yil 
1 19941 : lAbramowicj I1997I ) , flxing the viscosity parameter to 
Qv = 0.1, the ratio of gas to magnetic pressure to P — 10, 
and the fraction of turbulent energy in the plasma that heats 
the electrons to 5 — 0.3. The precise value of q for a partic- 
ular a ccreti ng black hole depends on m (see Equation 12 of 
lM03l ). lM03l show that, assuming the above parameters for an 
ADAF, sub-Eddington black holes with rh< 2 x 10"^ should 
have an average q = 2.3. Thus, in this work, we will use 
q = 2.3, dln^ B/dl nM = -0.5, and 91n <I>s/91n 771 = 0.5 
to compare to lMOSl 's theoretical FP coefficients for inverse 
Compton dominated X-rays. Differences in q between ob- 
jects (as well as the likely possibility that q varies with time) 
will contribute t o som e of the observed scatter about the FP. 

As noted in lMOSl . the above ADAF model represents a 
rather general description of a RIAF. Since the FP is a sta- 
tistical correlation, FP regression coefficients only provide 
information on the average properties of statistical samples 
of accreting black holes. In fact, every accreting black hole 
will have its own 'FP coefficients' that can vary with time, 
depending on its speciflc values of or, p, and q. Further- 
more, an individual black hole at a given accretion rate can 
even have different fractions of optically thin jet vs. inverse 
Compton emission contributing to the total observed X-rays 
depending on if it is h eading into outburs t or declining back 
into quiescence (e.g., iRussell et aLllioiol V We thus aim to 
keep our accretion model as general as possible so we can 
focus primarily on average trends. 

There are more complicated RIAF models in the liter- 
ature. For example, RIAFs can have convective instabilities 
if Ov becomes small enough (i.e., so-c alled convection dom - 
inated accretion flows or CDAFs; e.g, iNaravan et~alll2000l '). 
RIAFs can also produce strong mechanical outflows (i.e., 
advection do minated infl ow-outflow solutions or ADIOS; 
e.g., [Plandfor d fc Begelm an 1999). In such alternatives to 
ADAFs, X-ray emission is likely dominated by thermal 
bremsstrahlun^, which has q — 2 (and c|>r oc M~^7n, 
like a RIAF). For example, lYuan et al] l{20oi ) show that 
only for unphysical values of (5 = 1 (i.e., all the viscous 
dissipation goes toward heating the electrons) can ADIOS 
models predict that a radiative process other than thermal 
bremsstrahlung will dominate the observed X-ray spectrum 
of Sgr A* in the quiescent state. To capture the basic gen- 
eral characteristics of alternative RIAF models, we will also 
explicitly compare FP regressions to theoretical predictions 
assuming q — 2.0. We note though, according to Equation 12 
of lM03l . values of g ~ 2 could also be consist ent w ith the pure 
ADAF model described above (and used bv lMOSi ) if one does 
not consider any black holes with m < 10~*. Throughout 
this paper, we will therefore use the more general term RIAF 



We note, however, rapid X-ray variability (as is often observed) 
would usually rule out thermal bremsstrahlung. 
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Figure 1. Expected black hole mass (^a/) vs. radio luminosity 
Hr) coefficients from Equations [2] and [S] in the eases that X-rays 
are dominated by optically thin jet synchrotron (filled squares) or 
by inverse Compton off an X-ray corona (open symbols). For the 
corona model, we show predictions from a radiatively inefficient 
accretion flow (labelled as RIAF) with radiative efficieneies q = 2 
(triangles) and 2.3 (upside-down triangles). Each model is illus- 
trated by a pair of points connected with a dashed or solid line for 
the jet or RIAF models, respectively. These sets of lines/points 
show the range of expected coefficients for typical values of p, 
ranging from p = 2 (larger data points) to 2.8 (smaller data 
points). The top panel is for a radio jet with an inverted radio 
spectrum {an = —0.15), with the radio spectral index becoming 
increasingly steeper moving downwards. 



to describe inefficient accretion, and differences in RIAF 
models will be parameterized in a very general sense by the 
radiative efficiency q. 

The expected coefficients for jet synchrotron (Equa- 
tion [2]) vs. inverse Compton (Equation |3]) dominated X-rays 
are shown in Figure [1] (for different values of a_R, p, and g). 
Although different emission processes generally predict dif- 
ferent coefficients, the solutions are also mildly degenerate 
in that some combinations of models, ajj, p, and q can pre- 
dict similar coefficients. For illustrative purposes, we explic- 
itly show predicted FP coefficients for RIAF/corona models 
assuming q = 2.0 and 2.3. These two values represent (ap- 
proximate) extremes of a realistic range of plausible radia- 
tive efficiencies, as explained above. However, other values 
of q are also possible. 



3 FITTING THE PLANE: ARE FP SLOPES 
BIASED BY STATISTICAL EFFECTS? 

The observed multivariate correlation between multiwave- 
length emission and black hole mass supports the idea of 
black hole mass scaling, but the mere presence of a corre- 
lation on its own does not provide insight into the physics 
governing black hole accretion. Instead, one must compare 
the slope of the FP to predictions from physical models. 
Therefore, care must be taken when regressing the FP, to 
ensure the inferred slopes are not highly impacted by statis- 
tical effects. In this section, we perform a Bayesian regression 
analysis on a sample of accreting black holes, and we com- 
pare to the Merit function, which is the regression technique 
most commonly used by previous FP studies. 

3.1 Sam ple of Low-accretion Rate Black Holes 
from iKOel 

As stated earlier, low-accretion rate black holes across the 
mass scale are expected to be in a similar "hard-like" ac- 
cretion state, and they show the least scatter about the FP. 
Considering only low-accretion rate objects will therefore 
allow more sensitive constraints on how statistics may in- 
fluence the FP. Here, we use a sample of low- accre tion rate 
black holes assembled by lKOd which, following iKOfil . we refer 
to as the KFC sample|j 

The KFC sample includes radio luminosities, Lr = 
{i^Ljj)^Q^^, X-ray luminosities in the 0.5-10 keV band, and 
black hole masses for 77 accreting black holes. Included are 
25 observations of GBHs in the hard state (from GX 339-4, 
V404 Cyg, 4U 1543-47, XTE 1118-^480, and XTE J1550- 
564). There are 52 AGN, including one observation of Sgr A* 
('-lO'^ M0) during a hard X -ray flare, 17 low-luminosity 
AGN (LLAGN; ~10''-10** Mq), and 15 FR I galaxies and 
19 BL Lac objects (~10*-10® Mq). We use the KFC sam- 
ple because it contains exclusively low-accretion rate black 
holes, and its regr ession with various techniques has already 
been explored by IK06| . We exclude Seyfert galaxies, since 
their inclusion would increase the intrinsic scatter, and they 
may not be super massive analogs to ha rd state GBHs. We 
refer the reader to iPanessa et al.1 (|2007l ) for more informa- 
tion on_Seyferts regarding the FP. 

iKoel subdivide their sample into various subsets to test 
the effect of sample selection on the inferred FP slopes. We 
test the Bayesian regression on each of their subsets that 
includes GBHs and excludes Seyfert galaxies, leaving four 
KFC subsamples (see Table [ij which we refer to as the full 
sample (GBHs, Sgr A*, LLAGN, FR Is, and BL Lacs; 77 ob- 
jects), the contracted subsample (GBHs, Sgr A*, LLAGN; 43 
objects), the contracted+BL Lac subsample (GBHs, Sgr A*, 
LLAGN, BL Lacs; 62 objects), and the contracted+FR I 
subsample (GBHs, Sgr A*, LLAGN, FR Is; 58 objects). Al- 
though not discussed in this section, three other samples 
(called KFC-fSDSS, KFC-f SDSS-HBL, and KFC-fSDSS- 
LBL) are introduced in ^and included in Table [T] for com- 
pleteness. 

* The KFC sample is based on the sample compiled by IF04| . and 
it is modified to include more LLAGN and to exclude Seyfert 
galaxies. 
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Table 1. Samples Used to Regress the FP 



Sample 


Total 


Notes 




contracted 


43 


GBH, SgrA*, 


LLAGN 


contracted+BL Lac 


62 


GBH, SgrA*, 


LLAGN, BL Lac 


contracted+FR I 


58 


GBH, SgrA*, 


LLAGN, FR I 


full 


77 


GBH, SgrA*, 


LLAGN, FR I, BL Lac 


KFC+SDSS 


98 


GBH, SgrA*, 


LLAGN, SDSS BL Lac 


KFC+SDSS-HBL 


82 


GBH, SgrA*, 


LLAGN, SDSS HBLs 


KFC+SDSS-LBL 


59 


GBH, SgrA*, 


LLAGN, SDSS LBLs 



The KFC + SDSS, KFC + SDSS-HBL, and KFC + SDSS-LBL samples arc de- 
scribed in The other four are defined in il3. II 



3.1.1 Synchrotron Cooling 

The derivation of the FP in both IM03| and IF04| assumes 
synchrotron coohng is unimportant. The frequency where 
electrons undergo sig nifica nt c ooling is anti-correlated with 
black hole mass, and IF04| and IKOGI argue this frequency is 
below the X-ray band for the most massive black holes. Syn- 
chrotron cooling is thus a concern for FR I galaxies and 
BL Lac objects. Using the observed X-ray luminosities of 
FR I galaxies and BL Lac object is therefore not appropri- 
ate in the context of the jet model. Plus, X-ray emission may 
be SSC or external inverse Compton (EC) for some BL Lac 
objects (see ^4.3.2|l . The "X-ray luminosities" for the FR I 
galax;ies and BL Lac objects in the KFC sample are thus 
extrapolated from their optical nuclear luminosities, assum- 
ing an optically thin spectral index ax=0.6. We explore the 
impact of this assumption in Sections 13.4.11 and 14.41 Since 
synchrotron cooling is not important in the X-ray band for 
lower black hole masses, "real" observed X-ray luminosities 
are used for GBHs, Sgr A*, and LLAGN. 



3.1.2 Doppler Boosting 

Observed fiuxes of BL Lac objects are brightened by rel- 
ativistic beaming. Doppler boosting is negligible for other 
classes of AGN and GBHs considered here. To first order, 
the level of Doppler boosting for BL Lac objects is similar 
in the radio and optical bands (we extrapolate BL Lac X- 
ray fiuxes from the latter) , so beami ng do es not significantly 
impact the inferred FP coefficients ((FOJ). It may, however, 
contribute to so me but not all of the intrinsic scatter about 
the plane (e.g., iLi et al.ll2008l ). Since in this section we fo- 
cus on statistics, we do not correct for Doppler beaming 
yet. Such a discussion is deferred until ^4.4.11 We refer the 
reader to iHeinz fc Merlonil ()2004l l for a detailed discussion 
on how relativ i stic b eaming can affect the FP, as well as 
iGiroletti etHI lj2006l ). 



3.1.3 Error Budget and Intrinsic Scatter 

We adopt the error budget explained in detail in §2.3 of 
Black hole mass errors are typically of the order of 
-0.1 dex for GBHs, ~0.35 dex for LLAGN and FR I galax- 
ies (approximately the scatter about the Mbh — cr* rela- 
tion), and ~0.46 dex for BL Lac objects (which have more 
indirect mass estimates). Uncertainties on luminosities are 
dominated by errors in distance measurements, which typ- 
ically range from 0.1-0.4 dex. We also include <0.05 dex 



uncertainties on radio and X-ray fiux measurements. Addi- 
tional uncertainties contribute to the intrinsic scatter about 
the plane, ai^t ■ Potential sources of intrinsic scatter include 
the non-simultaneity of the AGN luminosity measurements, 
relativistic beaming, not all sources having identical accre- 
tion rates, the environments surrounding each black hole 
differing from source to source, and the radio and X-ray 
wavebands not probing identical regions of different source's 
SEDs. The last point means we do not observe exactly the 
same an and ax for every source, and we therefore do not 
expect every source to follow exactly the same FP coeffi- 
cients (see Figure [!}. Extrapolating "X-ray luminosities" 
from the optical for FR Is and BL Lacs abates this effect 
some, but not entirely (see tj4.4l for further discussion). 



3.2 Using the Merit Function to Regress the FP 

Fitting the FP is a multivariate regression problem, where 
a variable, j/i, depends on two independent variables, Xij, 
with j = {1,2}; each variable has a corresponding measure- 
ment error, ay. and cfx-^. Previous autho rs used the Merit 
function, a modified chi-square estimator (|Pressll2002l ): 



(4) 



The unknown parameters, Oj and 6, are found by minimiz- 
ing x'^- The c onsta nt b can be solved for analytically (see 
Equation 3 of IKOSI) . and the aj parameters can be found 
with a numerical optimization routine. 

Here we associate logLx with yi, logLn with xn, and 
logAf with Xi2, and the unknown parameters oi, 02 , and 
b correspond to ^r, ^a/, and B, respectively. IM03| assume 
measurement errors are isotropic (i.e., ctlx = — ""a/), 
and they adjust the uncertainties until = 1. This is an 
acceptable technique if the intrinsic scatter about the best- 
fit line dominates the error budget. However. IKOGI show this 
may be a strong assumption, and incorrect values for a, and 
b will be recovered if the error s are instead anisotropic. IK06| 
re-analyze the samples of lMOSi and lF04l . properly accounting 
for errors in each measured variable. They show the best-fit 
slopes of the FP strongly depend on the sample used, with 
the lowest-m AGN favoring the jet model. Including AGN 
with higher accretion rates increases the intrinsic scatter, 
and then the inferred slopes are perhaps more consistent 
with X-ray emission from inverse Compton. However, uncer- 
tainties on t he be st-fit coefficients from the Merit function 
regression m Hoi do not unambiguously identify the proper 
emission mechanism. 

The Merit function has some limitations. For example, 
it cannot handle coupled uncertainties, which is a potentially 
severe limitation since radio and X-ray luminosities are both 
dom inated by errors in distance measurements. Also, while 
iKOd properly account for measurement errors in each vari- 
able, their estimation of Oint is inferred by adjusting cri„t 
until is unity. We thus re-regress the FP with a more so- 
phisticated Bayesian technique (see next section) that can 
handle the above limitations (as wel l as l arge measurement 
uncertainties) to see if the results of IKOGI can be improved. 
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3.3 Bayesian Regression 

We use the method o f lKellvl (|2007l . hereafter IkOtI) to regress 
the FP. IKOTI take a Bayesian approach toward linear re- 
gression, estimating the probability distribution of the FP 
parameters given the measured datalf] The method assumes 
that the measurement errors are normal, that the intrinsic 
scatter about the FP is normal, and that the distribution of 
the independent variables can be approximated as a mixture 
of Gaussian functions. Under these assumptions it then uses 
a Markov chain Monte Carlo (MCMC) method to simulate 
random draws of the FP parameters from their probability 
distribution, given the measured data (i.e., the 'posterior' 
probability distribution). These random draws for the FP 
parameters then allow us to estimate their best-fit values, 
standard errors, and probability distributions. We refer the 
reader to IKOTI . especially their Sections 4 and 6, for more 
details. 

In this work, we typically take 10'' random draws from 
the posterior. Unless stated otherwise, we quote "best-fit" 
values of £,r, ^m, B, and aint as the median value of each 
quantity's posterior distribution. After checking that each 
posterior distribution follows an approximately normal dis- 
tribution, we report uncertainties in the text as the ilcr 
standard deviations . Since the best-fit regression coefficients 
are correlated, we will draw error ellipses in figures at both 
the Icr (68%) and 3a (99.7%) levels. 

In particular, this Bayesian method overcomes the two 
Merit function limitatio ns de scribed in the previous section. 
For one, the method of IKOTI can account for correlated er- 
rors. Secondly, it samples posterior distributions not just 
for the unknown parameters ^_r, £^m, and B, but also for 
the intrinsic scatter That is, aint is estimated directly 
from the data, and not by adjusting its value by hand until a 
modified chi-square estimator is close to unity. Especially for 
a correlation with relatively large intrinsic scatter like the 
FP, this has the attractive quality of not overweighting data 
points with small measurement errors that are h ighly scat- 
tered from the best-fit linelfl The method of lKOTi is also very 
well suited for handling heterogeneously selected datasets 
with large measurement uncertainties, like the KFC sample, 
as long as the intrinsic distributions of the independent vari- 
ables can be modeled as a mixture of Gaussians. This is a 
reasonable assumption as long as the number of Gaussians 
is large (the Gaussians do not need to be physically mean- 
ingful). Also, the measurement errors in the independent 
variables do not need to be of similar magnitude. 



3.4 Reliability of the Bayesian Regression 

Before direc tly c omparing Bayesian regressions usin g the 
technique of to the Merit function regressions in lK06l . 
we test the accurac y of the Bayesian technique on simu- 
lated FPs (also see IKOTI for more details on its reliability 
and comparison to other regression techniques) . We take the 
observed radio luminosities and black hole masses for each 



mliiunix.err .pro in the IDL Astronomy User's Library - 
|http : //Idlastro .gsf c .nasa. gov/ 1 

" Note, as for the Merit function, if measurement errors are over- 
or under-estimated, then cr^^^ is still under- or over-estimated, 
respectively. 




0.35 0.45 0.55 



Figure 2. Distributions of best-fit parameters with the Bayesian 
regression on 100 simulated FPs ( i|3.4l l. The input values (^/j, ^jv/i 
B, o-i„t)=(1.38,-0.81,-5. 00,0.40) are shown with dotted lines. The 
shaded regions illustrate the 16-84% confidence intervals {Icr). 



source in the full KFC sample (77 data points), assume these 
are the "true" values for each accreting black hole, and we 
predict "true" X-ray luminosities assuming {^r, ^m, B) = 
(1.38, -0.81, -5.00) and then add ±0.4 dex intrinsic scat- 
ter. The input FP slopes correspond to an = —0.15 and 
ax = 0.6 (i.e., p—2.2) in the jet model (see Figure [l|. Fake 
noise is then randomly added to each "true" value of Lr, 
Mbh, and Lx by assuming their errors are Gaussian with 
a standard deviation equal to the observed errors. When 
adding simulated noise, we take into account the fact that 
uncertainties in Lr and L x are correlated. We then regress 
the simulated FP with the lK07l Bayesian method, taking the 
best-fit parameters as the median values of the posterior dis- 
tributions for ^_R, S^M, B, and aint- We repeat this exercise 
100 times (i.e., we regress 100 simulated FPs), and the dis- 
tributions of the best-fit parameters are shown in Figure [5] 
The average best-fit parameters of the 100 simulated 
FPs match very well to the "true" input parameters, we 
recover {£_r) = 1.38 ± 0.04, (Ca/) = -0.81 ± 0.05, (B) = 
-5.12 ± 1.10, and {aint) = 0.42 ± 0.04 dex (median values 
of the fitted parameters are similar to the above means). 
The uncertainties quoted above are the standard deviations 
of the distributions shown in Figure [21 These uncertainties 
are consistent with the average of the la errors measured 
during the regression of each simulated FP: (ctj^ ) = ±0.04, 
((j^„) = ±0.06, (ffs) = ±1.16, and (cr<,.„,) = ±0.04. The 
Bayesian technique thus recovers accurate coefficients and 
intrinsic scatter about the FP, as well as realistic uncertain- 
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ties in each parameter provided our adopted error budget is 
reasonable. 



3.4-1 Synchrotron Cooling and Extrapolating FR I and 
BL Lac X-ray Luminosities From the Optical 

Extrapolating X-ray luminosities from optical nuclear lumi- 
nosities for FR I galaxies and BL Lac objects could system- 
atically bias the regressions. For example, the KFC sample 
assumes ax = 0.6, but using ax = 0.5 or q_r = 0.7 would in- 
crease or decrease the FR I and BL Lac "X-ray luminosities" 
by about a factor of two, respectively. To test the impor- 
tance of this effect, we take the 100 above simulated planes 
(which were simulated assuming an — —0.15 and ax = 0.6 
given the jet model). We then increase the FR I and BL Lac 
X-ray luminosities by a factor of two to simulate extrapo- 
lating their luminosities with ax = 0.5 and re-regress each 
simulated plane. We then repeat but decreasing their X-ray 
luminosities by a factor of two to simulate an extrapolation 
assuming ax = 0.7. The effect on the coefficients and 
^i\f from extrapolating with different values of ax is shown 
in Figure [S] Although the median best-fit coefficients of the 
100 simulations are formally consistent within iScr, larger 
values of ax appears to bias the FP toward shallower slopes. 
This is expected since steeper (i.e., larger ax) power laws 
will predict lower X-ray luminosities, so the correlation looks 
flatter. Similar intrinsic scatters, ai„t = 0.42-0.43 dex (with 
associated uncertainties of ±0.04 dex), are measured for all 
three values of ax ■ 

Finally, we simulate the situation where X-ray tele- 
scopes observe optically thin synchrotron jet emission with 
p = 2.2 (still assuming the radio is optically thick with 
an = —0.15) for lower mass (<10* -^q) black holes, but 
X-rays from higher mass black holes are synchrotron cooled 
and follow ax = 1- We test this by lowering the FR I and 
BL Lac X-ray luminosities by a factor of 12 and re-run 100 
simulations. From Figure [S] we see this pushes the coeffi- 
cients to even shallower slopes. We may have also measured 
an increase in intrinsic scatter, aint ~ 0.50 ±0.04 dex, which 
is not surprising. 

In ij4] we will confirm the assertions of IF04| and 
that extrapolating X-ray luminosities from lower frequen- 
cies is necessary if X-rays are dominated by optically thin 
jet synchrotron (i.e., one cannot test the theoretical FP co- 
efficients with observed X-ray fluxes for the most-massive 
SMBHs because synchrotron cooling is too strong). How- 
ever, our simulations show one must be careful in how this 
extrapolation is made. In the case of the KFC sample, if 
extrapolating with ax = 0.6, then the optical flux should 
be measured at a frequency where the non-thermal emis- 
sion actually follows a local power law spectrum ~ v~'^'^ , 
and therefore SED modeling is necessary (see We note, 
however, while this extrapolation effect is not negligible, in 
certain cases it may not be severe enough to strongly influ- 
ence the inferred X-ray radiation mechanism. For example, 
the KFC sample includes primarily low-accretion rate black 
holes with flat radio spectra an ~ —0.15. So, if regression 
of the KFC sample were to yield coefficients (and associated 
uncertainties) similar to our simulations, then it would likely 
still be possible to exclude the RIAF models in all four cases 
(especially the ones with very low-efficiency, i.e., q ~ 2.3). 
Regardless, this is a potential source of systematic bias one 
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Figure 3. The median value of the 100 best-fit coefficients in 
the Monte Carlo simulations assuming ax=0.5, 0.6, 0.7, and 1.0 
(cyan plus sign, blue filled circle, green diamond, and red cross, 
respectively) to extrapolate FR I and BL Lac X-ray luminosities 
from observed optical nuclear luminosities, la (solid blue line) 
and 3o" (dotted blue line) error ellipses are shown for ax = 0.6; 
errors are similar for the other three cases, but their error ellipses 
are not shown for clarity. Overplotted for reference are the the- 
oretical predictions for gii/ and ^/j for the optically thin jet syn- 
chrotron and RIAF models in the case aji = —0.15. The symbols 
showing the theoretical jet and RIAF predictions (i.e., squares 
and triangles connected by lines) have the same meaning as in Fig- 
ure[TJ. The simulated FPs assume a jet model with = —0.15 
and p = 2.2. Extrapolating X-ray luminosities from the optical 
can impact the best-fit coefficients. However, the effect (when 
Or = —0.15) is not severe enough to infer the incorrect radiation 
mechanism for X-rays from the best-fit coefficients, especially for 
very inefficient (i.e., large q) RIAFs. This figure appears in color 
in the online version of the article. 



should be concerned about, demanding special attention to 
how one analyzes multiwavelength data before inclusion in 
FP studies. 



3.5 Bayesian Regression of the KFC Sample 

Here, we apply the Bayesian regression on the contracted, 
contracted+BL Lac, contracted+FR I, and full KFC sam- 
ples (see tj3.1l and Table [T] for samp le definitions), and we 
compare to the regressions in IK06| using the Merit func- 
tion. Results are summarized in Table [2] and Figure |4] Save 
for the contracted+BL Lac sample, the Merit function re- 
gressions are consistent with the Bayesian regression at the 
3(7 level. However, there are important differences between 
the results from the two methods. First, the slopes recov- 
ered by the Merit function are consistently steeper than 
those from the IkotI method. We attribute some of this to 
the Bayesian regression properly handling correlated mea- 
surement errors between radio and X-ray luminosities. Also, 
the Bayesian regression recovers larger best-fit intrinsic scat- 
ters than the Merit function (note, the Merit function does 
not actually measure Oint directly from the data, rather its 
value is adjusted until ~ !)• Finally, uncertainties on 
the best-fit regression coefficients are consistently smaller 
with the Bayesian approach, which may eventually allow 
more definitive statements on the most-likely X-ray emis- 
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sion mechanism. We note, however, quoted uncertainties on 
the correlation coefHcients do not account for the potential 
biases from extrapolating FR I and BL Lac X-ray luminosi- 
ties from the optical discussed in ^3. 4. H and Figure [S] 

The best-fit regression coefficients from the Bayesian 
method for the contracted, contracted-f BL Lac and full 
KFC samples seem less sensitive to the subsample than for 
the Merit function. The regression coefficients for the con- 
tracted-f FR I subsample, however, are much shallower. The 
subsamples including FR I galaxies also have the largest in- 
trinsic scatter about the plane. We attribute this to several 
systematics related to difficulties in making the requisite ob- 
servations of FR I galaxies for FP studies (see H3.5.ip . That 
is, FR I galaxies bias the best-fit slopes (but the inclusion 
of BL Lac objects in the full sample seems to partially off- 
set this bias). FR I galaxies should be removed from the 
KFC sample until thes e systematics are more rigorously ad- 
dressed (although see iHardcastle et all |2009| regarding the 
placement of unbeamed radio galaxies onto the FP). 



3.5.1 Systematic Challenges to Including FR I Galaxies 

From a purely theoretical standpoint, FR I galaxies should 
follow the FP like other low-accretion rate black holes. They 
have sub-Eddington accretion rates, and, especially as m de- 
creases, their broad-band SEDs (i ncluding X-ray em ission) 
are generally jet-dominated (e.g., lEvans et al.ll200^ ). How- 
ever, there are several observationally based reasons to sus- 
pect that FR I galaxies could bias the FP regression. Most 
important is the effect of synchrotron cooling. For BL Lac 
objects, Doppler beaming will push the synchrotron cutoff 
toward higher observed frequencies, so that the optical wave- 
band still probes optically thin jet emission for many BL Lac 
objects. However, for unbeamed sources like FR I galaxies, 
even the optical could already be synchrotron cooled (also 
see H4.4p . Thus, the ax = 0.6 assumption for extrapolating 
FR I "X-ray luminosities" from the optical will introduce 
a larger bias than for BL Lac objects. Synchrotron cooling 
will thus systematically move the best-fit coefficients toward 
shallower values. Interestingly, the and ^a/ values of the 
contracted-f FR I subsample are similar to our Monte Carlo 
simulations where FR I and BL Lac "X-ray luminosities" 
are estimated using ax ~ 1 (to simulate synchrotron cooled 
jet emission in the optical. Figure [3]). 

Another potential source of systematic bias introduced 
by FR I galaxies is that it can be harder to isolate their 
nuclear radio emission, since their radio cores do not ap- 
pear as bright as BL Lac objects. The majority of FR I 

? alaxies in the KFC s ample have only VLA observations 
Chiaberge et al.lll999l '). meaning there is potentially some 
concern for contamination to their nuclear luminosities from 
optically thin radio emission that appears point-like at VLA 
resolution. (Of lesser, but perhaps not negligible, concern is 
radio emission from star formation contributing to the ob- 
served nuclear radio fiux). Thus, FR I radio spectra could 
be steeper than the rest of the objects in the KFC sample, 
which will also bias the best-fit regressions toward shallower 
slopes. Such systematics are likely less important than their 
synchrotron cutoff occurring below the optical waveband, 
especially since most of the KFC FR I galaxies are at low- 
redshift (2 < 0.1). We mention such potential systematics 



here primarily to illustrate that, in general, BL Lac objects 
provide "cleaner" probes of jet emission. 

Since BL Lac objects and FR I galaxies are presumed to 
be the same type of object but at different orientations, we 
prefer using BL Lac objects over FR I galaxies in this work. 
To properly include FR I galaxies, one would need to ensure 
observed radio core fiuxes are only including synchrotron 
jet emission, and SED modeling would be necessary to verify 
"X-ray luminosity" estimates are not infiuenced by dynamic 
cooling of the synchrotron emitting particles. 

3.6 Jet Synchrotron vs. Inverse Compton X-ray 
Emission 

The best-fit slopes of the FP can be used to investigate 
if X-ray emission from sub-Eddington black holes is dom- 
inated by optically thin jet synchrotron or inverse Comp- 
ton. However, we implicitly assume the jet interpretation 
when extrapolating FR I and BL Lac X-ray luminosities 
from the optical (such an extrapolation is improper if in- 
verse Compton). So we must omit the most massive black 
holes from this discussion, and we therefore only consider 
the contracted KFC subsample here. Even though the con- 
tracted subsample contains the smallest number of objects, 
and it has the most limited dynamic range in black hole 
mass and luminosity, the uncertainties on the Bayesian re- 
gression coefficients are narrow enough to make definitive 
statements on the dominant X-ray emission mechanism. We 
stress again that, in reality, both jet synchrotron and coro- 
nal inverse Compton components are likely present in all 
sub-Eddington black holes. The following discussion is thus 
aimed toward determining which emission mechanism, on 
average, dominates the X-rays observed from low-accretion 
rate black holes. 

Most (and perhaps all) accreting black holes in the KFC 
sample have an ^ 0, so we compare the contracted sub- 
sample regression coefficients to the predicted coefficients 
only for flat or inverted radio spectra. In Figure [5] we com- 
pare the best-fit regression coefficients to the coefficients 
that are predicted for aR — —0.15 and 0.00, in the cases 
where X-rays are dominated by jet synchrotron, or by in- 
verse Compton off a corona with radiative efficiencies q=2.0 
or 2.3. Considering even just la error bars on the Merit 
function regression coefficients, one cannot conclusively ex- 
clude any of the three models. However, the Bayesian re- 
gression favors the coefficient predictions if X-rays are dom- 
inated by optically thin jet synchrotron. As long as an ^0.0, 
the q=2.3 RIAF/corona model is excluded at the >3a (i.e., 
>99.7%) level. The Bayesian regression excludes the q=2.0 
RIAF/corona model with similarly high-confidence if the 
majority of objects in the KFC have inverted radio spec- 
tra (which is likely the case), but perhaps less significantly 
(p ~ 0.997) if the majority of objects have perfectly flat 
radio spectra (i.e., the Bayesian regression's 3a error el- 
lipse is similar to the expected RIAF coefficients if g = 2.0, 
an = 0.00, and p ~ 2.8 — 3.0). We thus conclude that the 
average black hole in the contracted KFC sample emits most 
of its X-rays as optically thin synchrotron radiation, and we 
exclude very radiatively inefficient RIAFs (i.e., q=2.3) with 
high- confidence. 

None of the theoretical predictions assuming aR — 0.5 
is consistent with the best-fit regressions, as expected since 
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Figure 4. Results of the Bayesian regressions (blue circles) compared to Merit function regressions (red crosses) from lKOa for the four 
KFC subsamplcs considered here. Icr (solid blue line) and 3o" (dotted blue line) error ellipses are drawn for the Bayesian regression, 
and la error bars are shown for the Merit function. Note, the Icr Bayesian uncertainties are smaller than the Ic Merit function 
uncertainties. The Merit function seems to systematically recover steeper slopes than the Bayesian regression. Synchrotron cooling moves 
the contracted+FR I best-fit slopes toward shallower values. This figure appears in color in the online version of this article. 



Table 2. Bayesian vs. Merit Function Regression of KFC Sample 



Sample 




Number 


Bayes 


an Regression 


Coefficients (tfiis 


work) 


Merit 


Function Coe 


fficionts [K06) 








B 


"int 




€a/ 


B 




Contracted 




43 


1.44 ± 0.09 


-0.89 ± 0.09 


-5.95 ± 2.58 


0.40 ± 0.06 


1.59 ± 0.21 


1.02 ± 0.21 


-10.15 ± 6.17 


0.12 


Contractcd- 


-BL Lac 


62 


1.42 ± 0.04 


-0.87 ± 0.05 


-5.37 ± 1.15 


0.34 ± 0.05 


1.64 ± 0.13 


1.08 ± 0.14 


-11.67 ± 3.59 


0.18 


Coiitractcd- 


-FR I 


58 


1.14 ± 0.06 


-0.60 ± 0.08 


2.65 ± 1.78 


0.48 ± 0.05 


1.25 ± 0.10 


0.74 ± 0.12 


-0.46 ± 2.93 


0.28 


Full 




77 


1.34 ± 0.06 


-0.81 ± 0.08 


-3.13 ± 1.58 


0.58 ± 0.05 


1.41 ± 0.11 


0.87 ± 0.14 


-5.01 ± 3.20 


0.38 



all KFC objects have flatter radio spectra. Also as ex- 
pected, radiatively efficient accretion flows are excluded. 
However, if we compare to predicted coefficient values as- 
suming an = 0.15, then the Bayesian regression analysis 
no longer excludes any of the three models. While the ma- 
jority of KFC black holes indeed have an < 0, a small 
number of LLAGN may have radio spectral indices extend- 
ing up to an = 0.3 based on non-simultaneous mult ifre- 
quency radio ima ging (see iTerashima fc Wilson! |2003| and 
iNagar et aLll2005l . from which the KFC LLAGN are taken). 
So, we searched NEeQ for all publicly available radio pho- 



tometric data for the 17 KFC LLAGN. We find only 4/17 
KFC LLAGN could potentially have an > 0, meaning >90% 
of the 43 objects in the contracted KFC subsample have 
Qij < 0. Thus, the presence of a small number of LLAGN 
potentially with or > may contribute to the intrinsic scat- 
ter about the FP, but they unlikely represent a large enough 
population to severely influence the best-fit (i.e., average) 
regression coefficients. We conclude that it is a reasonable 
assumption that our sample has flat/inverted radio spectra. 
However, it would be worth confirming our conclusion in fu- 
ture work with simultaneous multifrequency radio imaging. 



The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 



nology, under contract with the National Aeronautics and Space 
Administration. 
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Figure 5. Baycsian and Merit function regression coefficients for the contracted KFC subsample, compared to theoretical predictions for 
the jet (squares) and RIAF/corona models (triangles for q=2.0; upside down triangles for q=2.3) from Equations [2l and [3l We separate 
predicted FP coefficients for an = —0.15 and 0.00 into different panels for clarity. The Baycsian regression of the contracted subsample 
favors the jet synchrotron model. This figure appears in color in the online version of this article. 



It is important to once again stress that our conclusion 
favoring the synchrotron X-ray model applies specifically to 
the case of low-accretion rate black holes with flat /inverted 
radio jets. It is not appropriate to extrapolate this result to 
higher accretion rate SMBHs, which likely have X-rays dom- 
inated by emission from a corona. A similar study focusing 
only on black holes at higher accretion rates is out of the 
scope of this paper. 



4 FITTING VERY MASSIVE BLACK HOLES 
ONTO THE PLANE 

Our re-analysis of the KFC sample shows that one must 
overcome more systematics to obtain "X-ray" measurements 
from the most massive accreting SMBHs (>« 10* Mq). In 
this section, we shift our focus toward better understanding 
and quantifying these systematics. As described below, the 
inclusion of very massive SMBHs when regressing the FP 
will not yield additional insight into X-ray radiative pro- 
cesses from the lower mass black holes. However, many pre- 
vious FP studies include AGN with very massive central 
black holes, so an investigation focusing on potential biases 
introduced by those AGN is extremely important. 

We begin by removing the BL Lac objects from the 
KFC sample and we replace them with a uniform sam- 
ple of BL Lac objects fr om the Sloan Digital Sky Survey 

i SDSS; I^A et all [ 2000D with black hole mass estimates 
Plotkin et al]|201ll '). We then re-regress the FP with the 
same Bayesian technique as in the previous section. We ex- 
clude FR I galaxies in this section for the reasons described 
in H3.5.1I (also see H4.1|) . The Bayesian regression technique 
in conjunction with our recent and uniform SDSS BL Lac 
sample affords the requisite sensitivity to finally perform a 



rigorous study on biases introduced by the largest black hole 
mass bin. 



4.1 Why Use BL Lac Objects to Represent the 
High-end of the Mass Spectrum? 

First, it is important to reiterate the philosophical challenge 
to including very massive black holes in FP studies. In the 
previous section, we excluded the most massive KFC black 
holes (i.e., FR I galaxies and BL Lac objects) from our dis- 
cussion when attempting to determine the dominant X-ray 
radiation mechanism. The reason is that, for lower mass 
black holes, one can use real X-ray luminosities to differ- 
entiate between the optically thin synchrotron and inverse 
Compton scenarios. However, for the most massive black 
holes, synchrotron emission in the X-ray waveband is gen- 
erally radiatively cooled (or sometimes SSC/EC). So, even 
"jet" dominated SEDs from very massive black holes will not 
follow the theoretical predictions of Equation (2] if one uses 
"real X-ray" data. However, if very massive black hole SEDs 
are not jet dominated, then X-rays would still be sensitive 
to inverse Compton emission and real X-ray data should be 
used to place them on the FP. Thus, one must make an as- 
sumption on the dominant radiative process from the most 
massive black holes before deciding from which waveband to 
estimate their "X-ray luminosities." This a prion assump- 
tion on the radiative mechanism then makes it impossible to 
use the most massive black holes to diagnose the dominant 
radiation mechanism from a statistical black hole sample 
spanning the entire mass spectrumlf] 



Although theoretical FP coefficie nts for synchrotron cooled X- 
rays are predicted by iHeinj ||2004| '), comparing observations to 
those predictions presents a similar problem: observations in dif- 
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Even with the restricted dynamic range of the con- 
tracted KFC subsample, the Bayesian regression technique 
is sophisticated enough to infer that X-rays from those 
43 objects are generally dominated by optically thin syn- 
chrotron emission. To study in detail the systematics intro- 
duced by including the most massive black holes, we must 
add AGN with >^ 10* Mq central black holes back into 
the sample; we should only consider AGN that are very 
massive analogs to the types of accreting black holes in the 
contracted KFC subsample. BL Lac objects are the best 
very massive analogs : BL Lac objects hav e sub-Eddington 
accretion rates (e.g., iGhisellini et al.ll2009l ) and (uncontro- 
versially) jet-dominated SEDs with flat/inverted radio spec- 
tra. In principle we could also consider PR I galaxies, but 
we choose to restrict our massive SMBHs exclusively to 
BL Lac objects. The primary reason is the Doppler boosting 
of BL Lac jets: we can better isolate their jet core radio emis- 
sion (even without VLBI observations), and synchrotron 
cooling kicks in at slightly higher observed frequencies. The 
latter makes it easier to estimate optically thin synchrotron 
jet luminosities for BL Lac objects compared to PR I galax- 
ies. However, we refer the reader to lHardcastle et al.l (|2009l ) 
for a FP study that overcomes som e of the above sys t emat - 
ics: with X-ray spectral modeling, iHardcastle et all (|2009l ) 
are able to compare accretion-related vs. jet-related X-ray 
emission for a sample of unbeamed radio galaxies. 

In summary, by better controlling systematics at the 
high-mass end with a uniform sample of BL Lac objects, we 
can investigate how very massive black holes may bias FP re- 
gressions. Since the most massive SMBHs provide the widest 
dynamic, we will in the process also derive the most accurate 
FP coefficients to date in the context of the jet model. This 
assumption of the jet model is extremely important, as our 
handling of multiwavelength data for BL Lac objects is ap- 
plicable only if their SEDs are jet dominated. However, this 
assumption is also very reasonable in our case. We showed in 
the last section that the lower mass contracted KFC black 
holes are indeed dominated by synchrotron emission, and 
it is very well-established that BL Lac SEDs are jet dom- 
inated. That is, we know what BL Lac broadband SEDs 
ought to look like, allowing us to better investigate how the 
most massive SMBHs can bias FP regressions. 



4.2 Why Replace the KFC BL Lac Objects with 
Ones From the SDSS? 

Instead of using the BL Lac objects from the KFC sam- 
ple, we opt t o rep lace them with BL Lac objects from the 
SDSS. While iKOd do keep the KFC sample as uniform as 
possible, it is by nature heterogeneous, and this is especially 
true for their BL Lac objects. The 19 BL Lac o b jects in the 
KFC sample were taken from ISambruna et al. I (| 19961 ), and 
they were originally discovered either by the X-ray selected 
sample from the Einstein Observatory Extended Medium- 



Sensi tivity Survey fEMSS: iMorris et al.iri99ll : 



1991 ) or by the 1 Jy radio-selected sample (jStickel et al 



Stocke et al 



19911). These venerable BL Lac samples composed the first 



complete sets of BL Lac objects, and they provided tremen- 
dous advances toward our understanding of BL Lac phenom- 
ena. However, their shallow flux limits (and in the case of the 
EMSS, relatively small areal coverage) produced biased sam- 
ples that are not representative of the actua l BL Lac pop- 
ulatio n (e.g.. iLaurent-Muehleisen et al. l ll998l : |Plotkin" et al 



l2008h . Furthermore, only BL Lac objects for which F04 



ferent wavebands are necessary to uniformly probe synchrotron 
cooled emission across the entire mass scale. 



found black hole mass measurements in the literature made 
it into the KFC sample (and these black hole masses were 
also derived non-uniformly) . 

Since BL Lac objects provide such a large lever arm for 
regressing the FP, it is important to reduce as many of the 
above biases as possible to minimize the observational sys- 
tematics discussed in tj4.1l That is why we replace the KFC 
BL Lac objects with ones from the SDSS. Although BL Lac 
catalogs from the SDSS are not free of selection biases, these 
samples are large enough to include a popula tion of objects 
more represen tative of th e parent p opula tion (|Collinge et al.l 
[2005: Plotkin et al.ll201ol , hereafter [Piol ) . Perhaps more im- 
portantly in our opinion, the sample is uniform. That is, all 
objects are selected the same way, the multiwavelength data 
are taken with the same telescopes, and black hole measure- 
ments are derived with the same technique. 



4.3 The SDSS BL Lac Objects 

4.3.1 The KFC+SDSS Sample 

We use BL Lac objects from the l arge (7 23 ob ject) optically 
selected catalog from the SDSS (|PlO). IPIOI select BL Lac 
candidates based on their featureless optical spectra, apply- 
ing the standard criteria of Ca ll H/K breaks smaller than 
40% and no emission lines w ith rest-fram e equiv alent widths 
stron ger then 5 A (see, e.g.. lStocke et al..,1991. : iLandt et al.l 
|2002| ). For this FP study, we require radio and X-ray lumi- 
nosities, as well as central black hole mass measurements; 
we select a subset of 55 BL Lac objects for which we have 
data for all three FP axes. We refer to this sample as the 
KFC+SDSS sample (KFC GBHs, Sgr A*, LLAGN, and 
SDSS BL Lacs; 98 objects; see Table [TJ. 



4.3.2 The KFC+SDSS-HBL and KFC+SDSS-LBL 
Samples 

It is w ell known, e.g., from the blazar sequence (|Fossati et al.l 
Il998h . that BL Lac object SEDs can vary drastically from 
object to object, with their cutoff frequencies occurring 
anywhere from the near-infrared to the soft X-ray (e.g., 
iNieppola et al.ir2006l ). BL Lac objects with synchrotron cut- 
off frequencies in the near-infrared are referred to as low- 
energy cutoff BL Lac objects (LBLs) , those with soft X-ray 
cutoff frequencies are called high-energy cutoff BL Lac ob- 
jects (HBLs), and intermedia te-energy cutoff BL Lac obj ects 
(IBLs) are in between (e.g.. [Padovani fc Giomrnilll995l ). In 
addition to synchrotron emission, BL Lac objects also emit 
at higher frequencies via SSC and possibly EC processes. 
LBL/IBL X-ray emission is typically already SSC/EC (e.g., 
|Abdoelal][20lO), and it may be difficult (with extant data) 
for us to reliably estimate the luminosity of optically thin 
synchrotron emission for LBLs. However, this is much easier 
to do for HBLs, where X-ray emission is still synchrotron, 
albeit radiatively cooled. 
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With such a large parent sample of BL Lac objects, 
we can consider only HBLs and still retain a sufficiently 
large number of BL Lac objects. We thus also create the 
KFC+SDSS-HBL subsample, which consists of the 43 con- 
tracted KFC black holes and 39/55 SDSS BL Lac objects 
that we classify as HBLs (82 objects total; see Table [T|. We 
similarly create the KFC+SDSS-LBL sample, which only 
adds the 16 IBL/LBLs to the contracted KFC subsample 
(59 objects total; see Table[l|. We defi ne HBLs based on the 
ratio of X-ray to radio em is sion ( e.g., IPadovani fc Giommil 
1 19951 : iPerlman et"al] 1 19961 '). iPld followed relatively stan- 
dard convention, and they classified BL Lac objects with 
ctrx < 0.75 as HBLsH Here, we adopt a more conservative 
cut, and we identify our 39 HBLs as objects with a^x < 0.7 
(we calculate arx values from the rest-frame 1 keV X-ra y 
and 5 GHz radio luminosity densities in Table 7 of IPlOh . 
The remaining 16 BL Lac objects with arx > 0.7 are then 
IBLs and LBLs. The reason for our more stringent limit on 
X-ray brightness is to minimize the chance of an LBL or 
IBL contaminating the KFC-I-SDSS-HBL sample. Thus, we 
maintain with high confidence that the X-ray emission from 
these 39 HBLs is synchrotron emission and not SSC/EC. 
Comparing FP regressions of the KFC-fSDSS, KFC+SDSS- 
HBL, and KFC+SDSS-LBL samples will allow us to inves- 
tigate the impact of SSC/EC emission, as well as to fur- 
ther highlight the importance of comparing similar regions 
of SEDs across the entire accreting black hole mass spec- 
trum. 



4.3.3 BL Lac Black Hole Masses 

Quasar central black holes masses are normally estimated 
from the width s of broad emissio n lines and continuum lu- 
minosities (e.g.. lKaspi et aLll2000l ). However, the featureless 
nature of BL Lac objects that we exploit to isolate their non- 
thermal jet em ission also makes su ch measurements difficult 
(although see iDecarli et all I2OIII ). Instead, we limit our- 
selves to the low-redshift {z < 0.4) BL Lac subset that show 
a fiux component from their host galaxy in their SDSS spec- 
tra; the large size of our parent BL Lac sample ensures that 
a relatively large number of BL Lac objects show enough 
host galaxy flux. From the measured widths of their stellar 
absorption lines, we infer central black hole masses using the 
Mbh— o" relation (|Tremaine et al.ll2002l V We are able to mea- 
sure black hole masses for 71 SDSS BL Lac objects, with a 
precision ar ound 0.30-0.35 dex. These black hole masses are 
available in lPlotkin et al.l (|201ll ). where we also describe our 
technique in detail. Here, we include 55 of these 71 objects 
that also have X-ray detections in RASS, 39 of which we 
classify as HBLs (see ^4.3.2|) . We note that our requirement 
of substantial host galaxy flux biases our BL Lac subset with 
black hole masses toward the most weakly beamed objects, 
thus reducing (though not negating) the effect of Doppler 
boosting. Also, since all 55 of our SDSS BL Lac objects have 
z < 0.4, potential biases from luminosity evolution are mini- 
mized. This could be important because BL Lac objects may 



^ arx is the broadband X-ray to radio spectral index: arx = 
-log (^L^_i kfiv/L„,s GHz) /7.68, where L^_i ^eV and L^_5 qHz 
are X-ray and radio specific luminosities at rest-frames 1 keV and 
5 GHz, respectively. 



have a peculiar "negative" cosmic evolution (|Morris et al.l 



4.3.4 BL Lac Radio Luminosities 



iPld explicitly consider only optical properties for their 
BL Lac selectionl3 Post-selection, they correlated their 
sample to the Faint Images of the Ra dio Sky at Twenty 
cm survey (FIRST; iBecker et al]|l995l) and to the NRAO 
VLA Sky Survey (NVSS: ICondon et al.|[l99i ') radio surveys 
at 1.4 GHz , and to the ROSAT AU Sky Survey (RASS; 
IVoges et all 1 19991 . [jjOOf ) in the X-ray (0.1-2.4 keV). We 
use their published radio luminosities [vL^) at rest-frame 
5 GHz, which were estimated from FIRST/NVSS assuming 
a local spectral index of aj? = —0.27. All of the BL Lac 
objects considered in our FP work have flrm spectroscopic 
redshifts from host galaxy features for estimating luminosi- 
ties. We include in our error budget 5% uncertainties on the 
radio flux densities, and 5% uncertainties on distances (the 
latter is dominated by uncertainty in the Hubble flow). 



4.3.5 X-ray Luminosities 

Our goal is to better understand how the adopted method 
for estimating BL Lac X-ray luminosities may influence 
FP regressions. So, we estimate X-ray luminosities in three 
different ways to explore the effect. First, we take each 
BL Lac object's observed RASS position sensitive propor- 
tional counter (PSPC) count rate from 0.1-2.4 keV. From 
these count rates and each object's SDSS redshift, we es- 
timate broadband X-ray luminosities from 0.5-10 keV rest- 
frame, corrected for Galactic absorption, using the Portable , 
Interactive Multi-Mission Simulator (PIMMS : [Mukaill 19931 ) . 
Hydrogen colu mn densities along each object's sightline are 
taken from the IStark et all (1199^ 1 hydro gen maps using the 
colden tool in the Cha ndra Interactive Analy sis of Obser- 
vations package (CIAO: iFruscione et al.ll2006l ). We assume 
a local X-ray spectral index of qx = 1.50 and 1.05 for HBLs 
and LBLs, respect ively, which are typical va lues for each 
subclass (see, e.g., iPadovani fc Giommil 1 199^ ). We refer to 
these luminosities as "real X-ray luminosities." 

Uncertainties on the "real X-ray luminosities" are dom- 
inated by the adopted values of otx , and we conservatively 
add an uncertainty of ~25% to each "real X-ray luminos- 
ity" measurement. This error budget is adopted because 
that is the magnitude our "real X-ray luminosities" would 
change if we instead used local X-ray spectral indices be- 
tween 1 < ax < 3 (which is a reasonable range observed 
for BL Lac objects). We expect these "real X-ray luminosi- 
ties" to be affected by synchrotron cooling for HBLs, and 
dominated by SSC/EC for LBLs (and SSC/EC potentially 
contributing non-negligibly to the observed X-rays for IBLs). 
Thus, we do not expect FP regression coefficients using "real 
X-ray luminosities" to match the theoretical predictions in 
Equation [21 

We also estimate 0.5-10 keV rest-frame X-ray lumi- 
nosities from the observed optical flux densities assuming 
ax = 0.6 as in the KFC sample. We use optical luminosity 



However, multiwavelength constrai nts implicitly influe nce 
their selection for some sources (also see lCoUinge et al.|[2005l '). 
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densities of the AGN component (i.e., decomposed from the 
host galaxy) at rest-frame 5000 A. We refer to these X-ray 
luminosities as "KFC-like." The "KFC-like" luminosities are 
sensitive to the assumed value of ax = 0.6, since we are ex- 
trapolating these luminosities over 2-4 decades in frequency. 
We thus add a factor of two uncertainty to our "KFC-like" 
X-ray luminosities, the amount we expect luminosities to 
change if we rather choose ax = 0.5 or ax ~ 0.7 (uncer- 
tainties on distance measures are negligible compared to this 
factor of two uncertainty). Note, ax = 0.6 is a typical spec- 
tral index for optically thin synchrotron emission. So we do 
not expect this choice to systematically affect our best-fit FP 
(i.e., "average") regression coefficients; rather the fact that 
different objects have different ax values will primarily add 
intrinsic scatter. If the optical waveband really is probing op- 
tically thin synchrotron emission for BL Lac objects, then 
we expect the "KFC-like" luminosities to yield FP regression 
coefficients similar to the contracted KFC subsample. If the 
BL Lac synchrotron emission is already synchrotron cooled 
in the optical, then we expect shallower FP coefficients (see 

Finally, we build multiwavelength SEDs for each BL Lac 
object by correlating to other large-scale multiwavelength 
surveys, i ncluding the Westerb ork Northern Sky Survey 
fWENSS: [Rerigelink et"al]|l997l') a nd the Green Bank 6 cm 
survey (GB6: iGregorv et all Il996l }_ in the radio, the T wo- 
Micron All Sky Survey f2MASS: ISkrutskie et aLll2006l ) in 
the near-infrared, the Galaxy Evolution Explorer (GALEX) 
in the ultrav i olet, and the XMM-Newton Slew Survey 
jSaxton et all 120081 ') in the X -rays, in addition to the 
FIRST/NVSS radio and RASS X-ray data points avail- 
able from iPld . We fit parabolas (in log uF,j ~ log u) to 
these SEDs: loR'^J '^ = Ai{\ogiyf + Bilogv + Ci (e.g., 
iMassaro et al.ll2004l '). For the 39 HBLs we include X-ray data 
in the SED fit. For the other 16 BL Lac objects, it is unclear 
from their SEDs if the observed X-ray data are probing the 
synchrotron or the SSC/EC component, so we omit X-ray 
data from the fits for those 16 objects. 

From these SEDs, we measure the luminosity density 
at the frequency where the SED has a local spectral in- 
dex a = 0.6, corresponding to log = 0.4/(2y4i) + logVp, 
where logi'p = —Bi/2Ai is the peak frequency (i.e., where 
dlog vFv/d\og u = Q). log Up corresponds to the synchrotron 
cutoff frequency. We then extrapolate X-ray luminosities at 
rest-frame 0.5-10 keV from the flux density where a = 0.6. 
We refer to these X-ray luminosities as "SED-based." We 
add a factor of two uncertainty to the "SED-based luminos- 
ity" estimates. We expect the FP regression with "SED- 
based luminosities" to be consistent with the contracted 
KFC subsample regression. 

A potential advantage of the "SED-based luminosities" 
over the "KFC-like luminosities" is that the "SED-based lu- 
minosities" are always extrapolated from optically thin syn- 
chrotron emission, regardless of the frequency of the syn- 
chrotron cutoff. In principle, the "KFC-like luminosities" 
can probe different regions of each SED from object to ob- 
ject (i.e., some "KFC-like luminosities" could be synchrotron 
cooled), which would affect the best-fit regression coeffi- 
cients. We thus expect that using observed SEDs to estimate 
"X-ray luminosities" (even modeling with crude parabolas) 
will reduce the magnitude of this effect, provided our SED 
models are accurate parameterizations of each BL Lac SED. 



As noted at the beginning of this section, we restrict 
ourselves to 55 BL Lac objects with black hole mass mea- 
surements and X-ray detections in RASS (39 HBLs and 16 
IBL/LBLs for the KFC-^SDSS-HBL and KFC+SDSS-LBL 
samples, respectiv ely). In principle, we could potentially in- 
clude all 71 of the lPld BL Lac objects wit h bla ck hole mass 
measurements: the regression technique of IKOTI can account 
for upper limits when regressing the samples with "real X- 
ray luminosities." The other two types of "X-ray luminosi- 
ties" are extrapolated from lower frequencies, so their regres- 
sions do not require the handling of censored data. However, 
comparison between the three regressions with the three 
different "X-ray luminosities" are not as uniform if some 
include censored data and some do not. Perhaps more im- 
portantly, the X-ray data is very useful for constraining the 
HBL BL Lac SEDs and estimating accurate "SED-based X- 
ray luminosities" for HBLs. Reducing the above systematics 
provides a greater benefit than adding 16 more BL Lac ob- 
jects to the regression. Thus, we choose to restrict this study 
exclusively to SDSS BL Lac objects with RASS X-ray de- 
tections. 



4.4 Bayesian Regression Including SDSS BL Lac 
Objects 

The best-fit FP coeflicients for the KFC+SDSS, 
KFC-fSDSS-HBL, and KFC+SDSS-LBL samples are 
shown in Table [3 and the results for the KFC+SDSS-HBL 
sample are illustrated in Figure [S] Using real X-ray data 
clearly gives FP slopes shallower than any predictions from 
the optically thin synchrotron model. Thus, X-ray luminosi- 
ties indeed need to be extrapolated from lower frequencies 
for the most massive black holes, if their SEDs are jet 
dominated. The effect of SSC/EC emission contributing to 
the observed X-rays from LBLs is evident by comparing 
the best-fit "real X-ray" coefficients for the KFC-fSDSS, 
KFC+SDSS-HBL, and the KFC+SDSS-LBL samples. The 
inclusion of LBLs yields shallower coefficients. 

Interestingly, the best-fit "real X-ray" coefficients for 
the KFC+SDSS-HBL regression are nearly identical to the 
Bayesian regression coefficients for the contracted+FR I 
KFC subsample. Both samples therefore suffer from a sim- 
ilar type of bias due to the most massive black holes that 
is of comparable magnitude. Because we are certain that 
HBL X-rays are predominantly synchrotron cooled radia- 
tion (i.e., there is little to no contamination from SSC/EC 
or accretion fiow X-rays), emission in the optical waveband 
from the KFC FR I galaxies is also likely synchrotron cooled. 
Our FP simulations in i)3.4.1l further support the interpreta- 
tion that FR I galaxies emit synchrotron cooled radiation in 
the optical. Both sets of coefficients are consistent with the 
simulation where we modeled X-rays from lower mass black 
holes dominated by optically thin synchrotron {ax = 0.6) 
but X-rays from the most massive black holes (>~ 10* ^q) 
are synchrotron c ooled {ax = 1.0; see red cross in Figure [3]). 

iHeinj (|2004h show how to incorporate synchrotron cool- 
ing into the equations for scale-invariant jets by adding 
another scale length to the problem - the scaled-distance 
from the black hole where electrons become synchrotron 
cooled. The X-ray luminosity for a synchrotron cooled, scale- 
invariant jet should then follow Lx oc Mm, i.e., similar to 
that for a radiatively efficient {q = 1) standard accretion 
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Table 3. Baycsian Regression of the KFC+SDSS Samples 

Type of X-ray Beamed BL Lac Objects Debeamed BL Lac Objects 

Luminosity in B CT;„t Cj, Cjjj B CTi„t 

KFC + SDSS-HBL Sample (82 Objects) 

SBD-based 1.45 ± 0.04 -0.88 ± 0.06 -6.07 ± 1.10 0.07 ± 0.05 1.42 ± 0.05 -0.87 ± 0.07 -5.30 ± 1.52 0.07 ± 0.04 

KFC-like 1.43 ± 0.04 -0.88 ± 0.06 -5.71 ± 1.14 0.08 ± 0.04 1.41 ± 0.04 -0.86 ± 0.06 -4.90 ± 1.28 0.08 ± 0.04 

Real X-ray 1.16 ± 0.04 -0.62 ± 0.06 2.27 ± 1.10 0.45 ± 0.04 0.75 ± 0.07 -0.25 ± 0.08 13.85 ± 1.88 0.56 ± 0.05 

KFC + SDSS-LBL Sample (59 Objects) 

SED-basod 1.47 ± 0.04 -0.92 ± 0.06 -6.91 ± 1.19 0.16 ± 0.08 1.46 ± 0.05 -0.90 ± 0.07 -6.47 ± 1.50 0.16 ± 0.08 

KFC-like 1.35 ± 0.04 -0.81 ± 0.05 -3.24 ± 1.05 0.12 ± 0.07 1.32 ± 0.05 -0.79 ± 0.06 -2.30 ± 1.43 0.13 ± 0.07 

Real X-ray 0.96 ± 0.04 -0.43 ± 0.06 7.91 ± 1.23 0.52 ± 0.05 0.67 ± 0.07 -0.17 ± 0.09 16.13 ± 2.09 0.65 ± 0.07 

KFC-I-SDSS Sample (98 Objects) 

SED-basod 1.45 ± 0.04 ^-0.89 ± 0.06 -6.32 ± 1.13 0.11 ± 0.06 1.44 ± 0.05 -0.88 ± 0.07 -5.80 ± 1.45 0.11 ± 0.05 

KFC-like 1.41 ± 0.03 -0.86 ± 0.05 -4.92 ± 0.99 0.07 ± 0.04 1.37 ± 0.05 -0.82 ± 0.06 -3.77 ± 1.37 0.07 ± 0.04 

Real X-ray 1.06 ± 0.04 -0.51 ± 0.06 5.08 ± 1.14 0.55 ± 0.04 0.68 ± 0.06 -0.16 ± 0.08 16.04 ± 1.69 0.56 ± 0.04 
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Figure 6. Best-fit coefficients from the Bayesian regression of 
the KFC-f SDSS-HBL sample. For clarity, we do not show the 
best-fit coefficients for the KFC-I-SDSS or the KFC-I-SDSS-LBL 
samples (but see Table [S]) From top to bottom we show the 
best-fit regression coefficients when using SED-based (blue cir- 
cle), KFC-like (green circle), and real X-ray luminosities (red cir- 
cle) for the BL Lac objects, respectively, with Icr (solid line) and 
SfT (dotted line) error ellipses overdrawn. As in previous figures, 
the filled squares connected with a dashed line are the theoretical 
FP coefficient predictions for the jet model, for = — 0.f5 and 
2 < p < 2.8, as relevant to the KFC-I-SDSS samples. This figure 
appears in color in the online version of this article. 



disk[3 Cooling affects the dynamics of only the highest en- 
ergy relativistic ele ctrons, so th e expected radio emission 
remains unchanged. iHeinj \2004 ) then derive the expected 
FP coefficients for synchrotron cooled X-rays, and finds that 
one expects shallower slopes than if the X-rays are optically 
thin. Thus, our regression with "real X-rays" is biased to- 
ward the expected direction. However, it is not possible to 
direc tly compare our best-fit coefficients to the prediction 
from iHeinj |2004l ) more rigorously because the KFC X-ray 
luminosities for lower-mass black holes (i.e., GBHs, SgrA*, 
and LLAGN) are probing optically thin (i.e., uncooled) syn- 
chrotron in the context of the jet model. 

The "SED-based" regression coefficients for the 
KFC+SDSS, KFC+SDSS-HBL, and KFC+SDSS-LBL sam- 
ples are all consistent (within ilcr) with the regression coef- 
ficients for the contracted KFC subsample. We thus conclude 
that SED modeling allows accurate measurements for opti- 
cally thin synchrotron luminosities. We also find that the 
"KFC-like" regression of KFC+SDSS-HBL sample is con- 
sistent with the contracted KFC subsample. Thus, unlike 
for FR I galaxies, it is generally appropriate to extrapo- 
late HBL "X-ray luminosities" from the optical. That is, 
because BL Lac jet emission is Doppler boosted, the op- 
tical waveband typically probes optically thin synchrotron 
emission for HBLs. However, the "KFC-like" regression of 
the KFC+SDSS-LBL sample is not consistent with the con- 
tracted KFC subsample. We argue below that, like FR I 
galaxies, the discrepancy is because LBL synchrotron emis- 
sion (for at least some LBLs) is already synchrotron cooled 
in the optical waveband 

In Figure [7] we show the spectral index at rest-frame 
5000 A (i.e., where we measure optical luminosities for the 
"KFC-like" regression) inferred from our SED fits, and we 
compare to the frequency where each SED has = 0.6 
(i.e., from which we measure "SED-based" luminosities). All 
of the HBLs (filled circles) have spectral indices indicative 
of optically thin synchrotron emission (0.5 < < 1.0) 



Optically thin X-r ay jet emission would follow Lx oc 
(Mm)(i'^/i2)+(2/3)ax jMarkoflf et al.ll2003l : IfoI) . 
12 The "KFC-like" regression coeflScients for the KFC+SDSS 
sample naturally lie in between those for the KFC+SDSS-LIBL 
and KFC+SDSS-LBL samples (since the latter two samples are 
subsets of the KFC+SDSS sample). 
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Figure 7. The local spectral index of each SED fit at rest- 
frame 5000 A vs. log frequency where each SED fit has = 0.6. 
The filled circles are HBLs {oirx < 0.7), and the open circles are 
IBL/LBLs {oLrx > 0.7). The dotted vertical line marks 5000 A 
for reference. All HBLs have spectral indices at 5000 A indica- 
tive of optically thin synchrotron omission. However, synchrotron 
cooling may mildly bias the regression toward shallower best-fit 
coefficients if LBL/IBL X-ray luminosities are extrapolated from 
the optical. This effect is much more significant for FR I galaxies. 



at rest-frame 5000 Al3 Thus, the "KFC-like" and "SED- 
based" luminosities probe similar regions of HBL SEDs; so, 
we expect the FP regression coefficients of the KFC-f SDSS- 
HBL sample using those two luminosity estimates to then be 
similar and consistent with the contracted KFC subsample 
regression (and indeed they are almost indistinguishable). 

On the other hand, IBL/LBLs (open circles) have 
steeper spectral indices at rest-frame 5000 A, and many 
LBLs even appear to be radiatively cooled {a^, > 1). There- 
fore, the KFC-I-SDSS-LBL regression coefficients should be 
shallower than the KFC+SDSS-HBL coefficients. Thus, if 
one has the requisite data available to model SEDs, then 
that is a preferable method for estimating "X-ray lumi- 
nosities" from very massive accreting black holes with jet- 
dominated SEDs. Note, the "KFC-like" coefficients for the 
KFC-I-SDSS-LBL sample are not as shallow as the "real X- 
ray" coefficients for the KFC+SDSS-HBL sample because 
apparently not all IBL/LBLs are very strongly affected by 
synchrotron cooling at 5000 A rest-frame. 

From Figure [T] we can expand further on why FR I "X- 
ray luminosities" cannot simply be extrapolated from the 
optical. For unbeamed AGN with the most massive SMBHs 
(i.e., FR Is), the jet will appear to become optically thin at 
much lower frequencies (by almost an order of magnitude 
if BL Lac objects have Doppler pa rameters 5 ~ 7; also see 
Figure 6 of iBalmaverd e et al.ll2006l for a sketch of how FR I 
and BL Lac SEDs diflfer because of Doppler beaming). From 
Figure [T] we estimate a debeamed BL Lac object would 
have a steeper a,^ at 5000 A so that ai, > 0.8 always (and 
most with > 1.0). Thus, optical nuclear luminosities of 
FR Is should strongly be affected by synchrotron cooling. 



All of the HBLs also have synchrotron peak (i.e., cutoff) fre- 
quencies Up > lO^''" Hz, consiste nt with the expec ted cutoff fre- 
quencies for HBLs (e.g., see §6 of lAbdo et alll2010t) . 



SED-modeling of unbeamed jet-dominated AGN with very 
massive central black holes is thus necessary to place them 
onto the FP. 

The KFC-I-SDSS-HBL sample minimizes concern of 
synchrotron cooling systematically biasing the FP regres- 
sion. We thus consider the following regression to be the 
most robust: 

log Lx = (1.45 ± 0.04) log Lr - (0.88 ± 0.06) log Mbh 

-6.07 ±1.10, (5) 

To our knowledge, because of our sample selection and 
adopted regression technique. Equation [5] is the most ac- 
curate FP regression to date for sub-Eddington accreting 
black holes with flat/inverted radio spectra. 

For illustrative purposes, and comparison to previous 
FP studies, we show a projection of our final FP in Fig- 
ureH Shown is the best-fit for the KFC-f SDSS-HBL sample, 
with "SED-based" X-ray luminosities, both observed (top 
panel) and corrected for Doppler beaming (bottom panel; 
see H4.4.ip . For reference, we also show the location of FR I 
galaxies on the FP (with "X-ray luminosities" extrapolated 
from the optical), although they are not included in the fit. 
As expected, FR I galaxies tend to undershoot the FP. We 
note that regressing Lx and Ljj just for the SDSS BL Lac 
objects does not follow the same slope as the FP. This re- 
sult is due to the limited dynamic range when considering 
only the BL Lac objects, and the relatively large intrinsic 
scatter. Also, the dispersion in our SDSS BL Lac black hole 
mass measurements are of the same order as the measure- 
ment errors, ~0. 30-0. 35 dex. Thus, we do not expect to be 
able to infer reliable FP coefficients, especially for , when 
considering only BL Lac objects. 



4-4- 1 Doppler Boosting 

Throughout, we have assumed that Doppler beaming only 
affects BL Lac objects, and that, to first order, Doppler 
beaming does not largely affect the best-fit coefficients. That 
is, we assume that optically thick and optically thin syn- 
chrotron emission are beamed by approximately the same 
amount, so beaming primarily only moves BL Lac objects 
up or down the best-fit line in Figure[5](also see lFalcke et al.l 
2004'; 'Heinz fc Mcr lonil 12004 iLi et al.ll200i l. We note that 
Landt ct al. (2002) find no statistical evidence that more 
weakly beamed BL Lac objects (like the 55 BL Lac objects 
used in this work) have significantly different Doppler fac- 
tors in the radio and X-ray. Still, our assumption that the 
Doppler factor is similar in different wavebands is unlikely 
to be strictly true, and different Doppler factors probably 
contribute to some of the observed scatter. 

We test the effect of Doppler boosting by debeaming 
the SDSS BL Lac objects, assuming S — 7 \F04) . and thus 
reduci ng their radio and "X-ray luminosities" by a f actor 
5^+"" (|Lind fc Blandfordlll985l : IUrrv fc Padovanil(l995l ). We 
assume — —0.27 in the radio, Oi, = 0.6 for KFC-like and 
SED-based X-ray luminosities, and — 1.5 and 1.05 for 
HBL and LBL "real X-ray luminosities", respectively. The 
debeamed BL Lac objects have similar luminosities as FR I 
galaxies, as expected from the standard AGN unification 
paradigm (see Figure [Sjl. 

We then regress the KFC+SDSS, KFC+SDSS-HBL, 
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and KFC+SDS-LBL samples with the debeamed BL Lac lu- 
minosities. As long as our "X-ray luminosities" are probing 
optically thin synchrotron radiation, debeaming the BL Lac 
luminosities does not significantly change the best-fit re- 
gression coefiicients (see the coefficients for all three "SED- 
based" regressions and the "KFC-like" regression for the 
KFC-hSDSS-HBL sample in Table [S]). 

We thus conclude that Doppler beaming does not 
strongly infiuence the FP regression. Other systematics, 
such as making sure one probes consistent regions of ac- 
creting black hole SEDs across the mass scale, are likely 
more important. For example, when our "X-ray luminosi- 
ties" instead primarily probe synchrotron cooled emission 
or SSC/EC, as for the "real X-ray" regressions, then de- 
beaming the BL Lac objects does change the best-fit FP 
slopes. However, these regressions are not physically mean- 
ingful because we are probing different SED regions for dif- 
ferent black hole masses. Debeaming synchrotron cooled "X- 
ray luminosities" would move BL Lac objects down the syn- 
chrotron coo led radio/X-ray/mass correlation predicted by 
iHeinj (|2004t ). and not the "optically thin synchrotron" cor- 
relation followed by the lower-mass contracted KFC subsam- 
ple. Thus, debeamed synchrotron cooled or SSC/EC X-rays 
will bias the best-fit FP regression by a different amount 
than beamed X-rays. 



5 DISCUSSION 

5.1 Comparison of the Discovery Papers — Two 
Fundamental Planes? 

A few more comments on the two discovery papers and their 
infer red mechanisms for X-ray emission are in order. Because 
iMoal favors a RIAF and |FOJ a jet, there is sometimes an 
understandable but unfair perception in the liter ature that 
the two papers contradict each other. However, IM03| also 
outline how jet X-rays can lead to the FP, and in that case 
they predict coefficients similar to that of IF04| . There are 
likely at least two FPf^. and regressions of the FP only test 
the "average" emission mechanism. The situation is actually 
even more complicated, in that every source follows its own 
FP to some level due to different values of or, p, and q, 
which contributes to the observed intrinsic scatter about 
the FP. 

The apparently discrepant conclusions in the two dis- 
covery papers are an artifact of the samples used in each 
paper to test the FP. IF04| considered only low-accretion 
rate (i.e., jet-dominated) sourc es, while |M03| assembled a 
more diverse sample; the |M03| sample mixes some jet and 
some RIAF X-ray dominated black h oles, w hich naturally 
increases the intrinsic scatter. However, [mO^ exclude BL Lac 
objects, so their most-massive SMBHs include almost exclu- 
sively more luminous quasars that are well accepted to have 
coronal dominated X-ray emission. Given that those AGN 
provide the largest dynamic range and considering the pre- 
vious section's results, it is not surprising that the M03 re- 
gression favors the RIAF models (albeit with large intrinsic 

In fact, we have already alluded to a "third FP" in this paper 
from iHeinj | |2004) if synchrotron cooled X-rays arc considered 
instead of optically thin jet synchrotron. 



scatter) . That the FP coefficients are sensitive to the sample 
adopted in the two discovery papers (especially lower accre- 
tion rate objects seeming to favor jet X-rays), supports the 
idea that AGN mimic accre tion s tate s of GBHs. 

iKOel 's reanalysis of the lM03l andlF04l samples seems to 
confirm the conclusions o f the discovery papers. However, 
the coefficients found by IK06| for the two samples (using 
the same Merit function technique) are not rem arkabl y dif- 
ferent when considering the error bars; for the lM03l sam- 
ple, (CflLuna = 1-45 ± 0.17 and (^m)mo3 = -0-99 ± 0.22; 
for the IF04| sample, (5^)^04 = ^-^l ± 0.11 and i^M)„ r., = 
—0.87 ±0.14. It is thus an interesting "conspiracy" (as iKOg 
call it) that the two emission mechanisms yield FPs with so 
much overlap. This is in large part due to optically thin 
jet synchrotron emission and RIAF X-ray emission both 
approximately following L x ~ Mrri^ (|Falcke fc BiermannI 
1 19951 : iMarkoff et al.|[2003l ). Statistical techniques more so- 
phisticated than the Merit function are necessary to regress 
the FP with high enough accuracy to distinguish between 
the two radiation mechanisms. 

5.2 Comments on Other FP Studies 

5.2.1 I Yuan et al\ \200^) 

There is a well-established global radio/X-ray corre- 
lation followed by quiescent (i . e., hard state) GBHs : 
logLfl ~ 71ogL jf (|Corbel et al.ll2003l : ICallo et al.ll2003l ). 
lYuan et al.l (|2009r ) claim this correlation is best explained if 
the X-rays are pr edominantly from th e accretion flow. This 
claim is based on lYuan fc Cuil (|2005l ). who model emission 
from accreting black holes with components from both a 
jet and from an ADAF. Their ADAF accounts for the ef- 
fects of outflows and convective instabilities, and it includes 
emission from synchrotron and br emsstrahlung proce sses 
and their thermal Comptonization. lYuan fc Cuil (|20Q5l ) ar- 
gue that only when the X-ray luminosity drops below a crit- 
ical X-ray luminosity, Lx,crit ~ 10~^- 10~^ LEdd, can jet 
synchrotron instead dominate the observed X-ray emission. 
They predict that the logLa --^ 0.7 log correlation will 
then steepen to logLij ~ 1.23 log Lx, as the X-ray emission 
swit ches from being A DAF to jet dominated 

lYuan et al.l (|2009l ) use a sample of 22 very low-accretion 
rate AGN with Lx < 10~® Lsdd (i.e., be low Lx.cru) 
to reg ress the FP and test the prediction of lYuan fc Cuil 
(|200a ). They find logLfl ~ (1.22 ± 0.02) log Lx + (0.23 ± 
0.03) log Mbh, in excellent agreement with lYuan fc Cuil 
(|2005h R lYuan et al.1 (|2009l ) do not include any AGN with 
Lx ^ Lx,crit, but they note that the accreting black 

No te, di fferent from t he co nvention in this paper. pYuan &: Cuil 
||2005|) and lYuan et al. cast Ln as their dependent vari- 

able. So, their prediction for a steeper correlation at lower X-ray 
luminosities corresponds to a shallower (i.e., lower) value of ^/j in 
our notation using log Ly as the dependent variable. 

It is surprising that lYuan &: Cuil ||2005| ) report u ncerta inties 
on their coefficients smaller than both our work and |M03|. smce 
the latter two studios adopt much larger samples (plus wo use 
an arguably more acc urate regression method). Comparing to 
iCultekin etal] ll2009al) . who use d a similar r egress io n technique 
on a similarly sized sample as I Yuan et al.1 | |2009|) . I Yuan et al.l 
seem to underestimate errors on their coefiicients by an 
order of magnitude. 



© 2011 RAS, MNRAS 000.[Tlf23l 



18 Richard M. Plotkin et al. 




Figure 8. Our best-fit FP for low-accretion rate black holes (the KFC-I-SDSS-HBL sample; 82 objects) using the Bayesian regression 
algorithm and SED-based X-ray luminosities. The top panel shows the regression for beamed BL Lac objects, and the BL Lac objects 
are debeamed in the bottom panel. FR I galaxies are shown for reference, but they are not included in the regression. This figure appears 
in color in the online version of this article. 



holes in the sample used by lM03l all have Lx > Lx,crit- 
find a shallower FP regression for their more lumi- 
nous accreting black holes [logLij ~ (0.60to';[i) log Lx + 
(0.78lg jg) log Mbh], which is consistent with the radio/X- 
ray correl ation for quiescent GBHs (i.e., l ogLj; ~ 0.7 l og Lx ; 
iGallo et al. l()Oa) . The conclusion of lYuan et all (|2009l ') 
therefore supports lYuan fc Cuil (|2005l ). that the radio/X-ray 
correl ation of quiescent black holes described in lGallo et al.l 
is due to the jet in the radio and inverse Compton in 
the X-ray. However, below a critical luminosity, the radio/X- 
ray correlation steepens as the X-ray emission instead be- 
comes dominated by jet synchrotron. 

The argument in lYuan fc Cuil (|2005l ) for a steeper corre- 
lation when Lx < Lx.crit assumes that jet X-ray emission 



scales linearly with M. Then, Lx.crit corresponds to the 
critical accretion rate where jet X-rays start to outshine the 
accretion flow (RIAF X-ray emission scales approximately 
quadaratically with A'l). Their assumed jet X-ray depen- 
dence, Lx PC M, implie s that the jet X-rays are synchrotron 
cooled (see lHeinj[2003 '). On the other hand, if jet emission 
is still optically t hin in the X-rays for GB Hs, then Lx ~ 
^(i7/i2)+(2/3)(.x (|Markoff et al.ll2003l : [Foi ). Then, optically 
thin synchrotron can dominate over coronal em ission in the 
X-ray s at higher values of Lx/LEdd, and the iGallo et al.l 
(|2003l ) correlation can be interpreted as predominantly op- 
tically thin jet synchrotron X-rays. There is evidence that 
the X-ray flux of the GBH XTE J1550-564 is 100% optically 
thin jet synchrotron at Lx ~ 4 x 10^^ — 4 x 10^* LEdd as 
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it fades back into t he low-hard state following an outburst 
l|RusseU et al.ll20ld )F1 This X-ray luminosity is higher, al- 
though n ot totally in consistent, with the value of Lx,crit pre- 
dicted bv lYuan fc Cui (200^). H owever, an important differ- 
ence IS that iRusseli et all |2010l ) measure an X-ray spectral 
index ax ~ 0.7, indicating optically thin (i.e., uncooled) X- 
ray emission accounts for 100% of the X-ray flux of XTE 
J1550-564 at energies of a few keV (i.e. X-ray luminosity 
should scale approximately quadratically with M and not 
linearly, inconsistent with the assum ption of Yuan fc Cuil 
[iooi). It is also important to note that iRussell et al.l (|2010l ) 
find synchrotron emission cannot contribute all the hard X- 
rays at the beginning of the outburst, so at least one other 
parameter besides accretion rate is likely important for con- 
tro Uing the X-ray emission mechanism. 

iRussell et al.l (|2010l ) note that the synchrotron cutoff 
is believed to occur at >40-100 keV for hard-state GBHs, 
above the observed hard X-ray band (see their §4.1). We 
might t hen expect, if the quiescent GBHs used to mea- 
sure the lCallo et al.l (|2003l ) radio/X-ray correlation were ob- 
served at X-ray energies above the cutoff, then they would 
follow a correlation closer to the log Lu ~ 1.23 log Lx corre- 
lation predicted bv lYuan fc Cuil (|2005l ). That is, the change 
in slope does not necessarily require a switch between RIAF 
and jet dominated X-rays. Rather, it could instead more 
heavily depend on the waveband in which "X-ray luminosi- 
ties" are measured. 

We argue that it is this type of observation al effect re- 
spons ible for the diff erent FP slopes measured bv lYuan et al.l 
l|2009t ) compared to IM03| . More pre cisely, the X - ray ob ser- 
vations of the AGN sample used by lYuan et al.l (|2009l ) are 
affected by synchro tron cooling, while very few AGN in 
the sample used by IM03| are emitting pre dominantly syn- 
chrotron cooled X-rays. IVuan et al.l (|2009l ') use the lowest- 
luminosity AGN yet in any FP study, but most have such 
low Lx/LEdd because they have relatively large black hole 
masses. Of their 22 AGN, 14 have M > 10* Mq. Thus, ff 
their X-rays are indeed due to jet emission, they will be syn- 
chrotron cooled, an d lower best-fit val ues of (as predicted 
by lHeinj|2004l and lYuan fc Cuil l2005l) are ex pected. 

The premise behind lYuan et all (|2009l ) is valid, that 
synchrotron cooled X-ray emission can outshine the corona 
if Lx/Lsdd is low-enough, and indeed their 22 AGN seem 
to have X-rays dominated by cooled sync hrotro n emission. 
However, comparing their results to the lM03l coefficients 
is not evidence that all accreting black holes with L x > 
Lx,crit have coronal dominated X-ray emission. lM03l omit 
Lx > Lx,crit AGN with jet-dominated SEDs (i.e., BL Lac 
objects) from their sample, and this infiuences their regres- 
sion toward FP coefficients predicted by the RIAF model 
used by lM03l . Based on our results in EH ff lMOal had in- 
cluded BL Lac objects (with Lx > Lx.cri t) and used "real 
X-ray luminosities" for them, then lM03l also would have 
recovered coeffici ents more simila r to those recovered by 
lYuan et al.l (|2009l ') (e.g.. Figure O. IM03| note their FP scal- 
ings are not valid if synchrotron cooling becomes important, 
so one must be very cautious comparing the low-luminosity 



AGN in lYuan et all (|2009l ) to lMOl as such a comparison is 
not uniform. 

Ide Gasperin et al] (|201ll ) report on sensitive radio and 
X-ray observations of 16 Type 2 LLAGN from the SDSSH 
Their LLAGN sample covers a large dynamic range in lu- 
minosity, but the sample is selected to cover a narrow range 
in black hole mass (10** < Mbh < 10** '' Mq). From their 
radio, optical, and X-ray data, they conclude that these 
LLAGN have radio and X-ray em ission dominated by jet 
radiation. Ide Gasperin et al.l (|201ll ) also use their sample to 
regress the FP, and they find that the FP coeffici e nts o f 
their sample are consistent with th at o f lYuan et al.l (|2009l ) 
(and inconsistent with both IM03| and K06 ). Howe v er, al l 
of the LLAGN in the sample of Ide Gasperin et al.l (|201ll ) 
have X-ray luminosities brighter (most by more than an or- 
der of magn i tude) than the value of Lx,crit predicted by 
lYuan fc Cuil (|2005l ') . These higher lumino sities make i t chal - 
lenging to interpret the consistency with lYuan et al.l (|2009l ) 
being due to a transition from corona to jet dominated X-ray 
emission below a critical X-ray luminosity. Again, we argue 
tha t such an interpretation i s not necessary. The LLAGN in 
the lde Gasperin et al.l (|201ll ) sample have relatively massive 
central black holes, and their X-ray emission is very likely 
affected by synchrotron cooling. Thus, the FP slopes recov- 
ered bv lde Gasperin et al.l (|201ll ) are naturally explained if, 
because of their very massive central black holes, their X- 
ray observations (from 2-10 keV) are probing primarily syn- 
chrotron cooled jet X-ray emission (opposed to optically thin 
synchrotron that would be probed by X-ray observations if 
the AGN had lower cen tral black hole masses). 

iGallo et al.l (|2006l ) obtained simultaneous radio and 
X-ray observations for the GBH A0620-00 in quiescence 
with extremely low Lx ~ 10~*'^ Lsdd- They show that 
the observed radio luminosity of A0620-00 is consistent 
with extrapolating the logLij ~ 0.7 log Lx correlation to 
low Lx , inconsistent with the lYuan fc Cuil (|2005l ) predic- 
tion that the correlation should steepen below Lx,crit = 
lO"'' - 10"** Lsdd- It is true that the GBH radio-X-ray 
luminosity correlation has some intrinsic scatter, and one 
cannot definitively extend to very low accretion rates based 
on only one data point. However, if the correlation steep- 
ens to Lr ~ L3^2^ at Lx,crit < 10-^-^ LEdd, then A0620-00 
should be radio-fainter than the \og Lr ~ 0.7 log corre- 
lation by a factor of 39. This is larger than even a generous 
o"int = 1.00 dex scatter in radio luminosity. 

Finally, FR I and II radio galaxies may yield additional 
insight. X-ray emission in FR I galaxies is likely dominated 
by the jet, while more luminous FR II galaxies primarily 
emit coronal X-rays (see ii3.5.1l and references therein). The 
so-called FR I/II dichotomy (e.g., iLedlow fc Owed Il996l ) 
could be due to a state transition analog ous to that of 
hard and soft-state GBHs; IWold et all (120071 1 argue the state 
transition to jet dominated X-rays occurs near m = 0.004 
and is a factor of 10 lower for FR I/IIs than for GBHs. 
Estimates of a state transition in blazars (i.e., between 
BL Lac objects and fiat spectrum r adio quasars) are sim - 
ilarly placed around Q.QlLsdd (e.g.. [Ghisellini et al.|[2009l ). 



Their sample actually includes a 17*'' AGN, but they note that 

^"^ We convert the bolometric luminosities quoted in lRussell et al.l that source is most likely a BL Lac object and we thus exclude it 

||2010| ') to that in the 2-10 keV X-ray band. from this discussion. 
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Scaling downward, then we speculate it is reasonable to 
think sub-Eddington GBHs could also show jet dominated 
SEDs below a few percent LEdd- While further study is 
needed, we conclude based on all the evidence presented in 
this subsection, it is very likely that some accreting stellar 
mass black holes with Lx > W~'' — 10~^ Lsdd can have jet 
dominated SEDs in the X-ray waveband. 



emitting primiarly synchrotron cooled jet emission in the X- 
rays, while the RIAF is predominantly responsible for the 
radio-quiet objects' X-ray emission. First, however, further 
study into potent ial selection biases is necessary. For ex- 
ample, due to the [Anderson et al.l |2003) AGN being X-ray 
selected, perhaps the requirement of X-ray emission from 
RASS preferentially excludes radio-loud AGN at low-radio 
flux densities. 



5.2.2 \Li et al\ \200^) 



iLi et al.1 (|2008h examined the FP by regressing a subset 
of X-ray selected broad-line AGN (i.e., excluding BL Lac 
objects) from the SPSS and R ASS (their AGN are taken 
from the [Anderson et al.l [2OO7I X-ray AGN catalog) Of 
th e 10^ spectroscopically confirmed X-ra y emitting AGN 
in [Anderson et al.l (|2007r i. iLi et al.l (|2008l 'l use 725 in their 
study for which they can measure black hole masses from 
widths of broad emission lines. Their sample is uniform, 
and to our knowledge, the larg est number of o bjects imple- 
mented in a FP study to date. [lT et al.1 (|2008D find negligi- 
ble dependence on black hole mass, which they attribute to 
using different rest-frame radio and X-ray luminosities than 
IM03| (as well as to the heterogen eous n ature of the black hole 



mass measurements adopted by MOsI) . A pe rhaps more im- 
portant consideration is thatllT et al.1 (|2008D only use SDSS 
AGN with large black hole masses in their regression (each 
with relatively large uncertainties). There is thus unlikely a 
large enough dynamic range in black hole mass to reliably 
recover ^Af. From their Figure 1, their black hole masses 
are p eaked near 10*-^-10^ ° Mq; also see lKellv fc Bechtoldl 
I2OO7I on estimating AGN central black hole masses from 
single epoch spectroscopy of broad emission lines, and how 
the inferred distribution of mass measurements may appear 
bro ader than in real i ty. 

lAnderson et al.1 (|2007l 'l also include ~250 BL Lac ob- 
jects in their AGN catalog, which are excluded from the 
iLi et all (|2008l ') subsample since they do not display strong 
enough emission lines to measure black h ole m asses. The 
BL Lac selection algorithms implemented in lPlOl . our parent 
sample of BL Lac o bjects in this work, wer e largely based on 
those developed in [ Anderson et al.l l|2007l l (combined with 
Ande rson ct al.l l2003h ." Thus, the unbeamed AGN used by 
Li et al. (20ol) cover similar dynamic ranges in luminosity 



and black hole mass, with similarly sized error bars, as the 
beamed AGN used in our work. We find it is not possible to 
obtain reliable coefficients considering only SDSS BL Lac ob- 
jects. Rather, massive accreting black holes should be used 

in conjunction with their lower-mass counterparts. 

Finally, we note an interesting result that iLi et al.l 

l|2008l ) find different correlations for radio-loud and radio- 
quiet SDSS AGN [Lr oc ijf* if radio-loud, and Lr oc L^x' if 
radio-quiet). Given the large black hole masses of the sam- 
ple used bv iLi et all (|2008l '). their radio-loud AGN X-ray 
luminosities (as observed by RASS) would be affected by 
synchrotron cooling if jet dominated. The different corre- 
lations may then be consistent with their radio-loud AGN 



13 iLi et al.l 1 I2OO8I ) is an updated version of lWang etHI ||2006|") to 
include a larger X-ray selected AGN sample. IWang et al. I 1I2OO6I ') 
used AGN fr om a similar, bu t smal ler, X-ray selected SDSS AGN 
sample from lAnderson et al. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we explore how statistical effects and sample 
selection affect the best-fit slope of the FP. Since different 
coefficients are expected depending on the physical condi- 
tions and geometry very close to the black hole, and whether 
X-rays are dominated by optically thin jet synchrotron or 
inverse Compton, the FP can be exploited as a tool to di- 
agnose radiative processes of accreting black holes. Previ- 
ous FP studies have been very insightful. However, by us- 
ing a more sophisticated (Bayesian) regression technique, it 
is possible to refine the FP even further. Uncertainties on 
our best-fit coefficients are generally smaller than other FP 
studies, allowing more decisive statements on radiative pro- 
cesses. Also, our technique provides the first estimates of 
the intrinsic scatter about the FP measured directly from 
the data. 

Ideally, one desires the largest possible dynamic range 
in mass and luminosity to obtain the most reliable FP coef- 
ficients. We discuss how including only one class of AGN in 
FP studies severely limits the dynamic range and potentially 
decreases the reliability o f the infe rred c oefficients. However, 
we confirm the claims of IF04| and lKod that jet X-rays from 
the most massive SMBHs (i.e., >~10* -^o) ^^e strongly af- 
fected by synchrotron cooling. This fact makes it difficult to 
incorporate the most massive SMBHs in FP studies and then 
uniformly test if X-ray emission is dominated by the jet or 
by the corona. For the most massive black holes, real X-ray 
data should be used to test corona models; however, if jet 
dominated, then "X-ray luminosities" must be extrapolated 
from lower-frequency wavebands (to be sure one is compar- 
ing only optically thin jet synchrotron emission across the 
entire mass scale). 

An advantage of our Bayesian regression analysis is that 
it is robust enough to measure reliable FP coefficients even 
if the most massive SMBHs are excluded. Thus, we can 
test the dominate X-ray radiation mechanism without mak- 
ing any a priori assumption on the radiation mechanism. 
We find that a subset of sub-Eddington black holes with 
Mbh < 10* Mq (i.e., the contracted KFC subsample) favors 
optically thin jet synchrotron X-rays, provided their radio 
spectra are flat/inverted. Our analysis excludes with high 
probability (p > 0.997) the possibility that sub-Eddington 
black holes emit X-rays from a corona surrounding a RIAF 
with very low radiative efficiency (i.e., q = 2.3), or surround- 
ing a radiatively efficient q — 1 accretion flow. 

We add a subset of very massive accreting black holes 
with uncontroversially jet dominated SEDs (i.e., BL Lac ob- 
jects) to our sample of M < 10* Mq sub-Eddington black 
holes. For the reasons described above, the inclusion of these 
BL Lac objects does not provide new constraints on the 
X-ray radiation mechanism. Rather, we use these BL Lac 
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objects to illustrate the systematic challenges to including 
the most massive black holes in FP regressions. In partic- 
ular, we show how the effects of synchrotron cooling from 
jet dominated SEDs are increasingly important in the X- 
ray waveband (and even in the optical waveband for FR I 
galaxies and some LBLs) as black hole mass increases; we 
show that synchrotron cooling will bias FP regressions if not 
taken into account. We argue that one should model SEDs 
to measure the luminosity of their optically thin jet emis- 
sion, although for HBLs one can alternatively extrapolate 
optically thin synchrotron luminosities from the optical. 

The inclusion of BL Lac objects increases the dynamic 
range in black hole mass and luminosity, allowing more ac- 
curate estimates of the FP coefficients especially when we 
only add HBLs to the lower-mass black hole sample. For 
example, the uncertainties on our KFC-I-SDSS-HBL regres- 
sion coefficients are cr^^ = ±0.04, a^^^^ — ±0.06, compared 
to cr^fl = ±0.09, a^^, = ±0.09 for the contracted KFC 
sample. We find logi^ = (1.45 ± 0.04) log Lh - (0.88 ± 
0.06) log Mbh-6.07± 1.10, with cri„t = 0.07±0.05 dex, and 
a similar regression is obtained if we debeam our BL Lac 
objects. 

Our uniform SED modeling of a large number of BL Lac 
objects, combined with our adopted Bayesian regression 
technique, makes this one of the most accurate FP regres- 
sions yet performed. Note, the above correlation cannot be 
applied to higher accretion rate black holes, or to black holes 
without fiat or inverted radio spectra. These black holes may 
instead emit coronal X-rays and follow a different relation. 
We also argue that, because of several observational effects 
due to relativistic beaming, inclusion of BL Lac objects al- 
low more reliable coefficients than including FR I galaxies, 
which may also have jet dominated SEDs. On purely obser- 
vational grounds, we advocate the removal of FR I galaxies 
from FP studies unless these difficulties are rigorously ad- 
dressed. 

A challenge to fitting the FP is that every object fol- 
lows its own FP in principle (i.e., depending on its particular 
values of an, p, and q). Also, accreting black holes are vari- 
able: their accretion rates, or, p, and q values, etc., change 
with time, meaning even individual sources are expected to 
follow different FP correlations over time. Thus, statistical 
regressions only teach us about the "average" black hole. 
Complicating the matter even more is that a magnetized 
corona may not be a totally distinct entity fro m a coUimated 
magn etized outfiow at the base of a jet (e.g.. [Markoff et al.l 
[2003), and observational signatures differentiating the two 
mechanisms are rather subtle. In reality, both coronal and 
optically thin jet emission likely contribute to the observed 
X-ray emission, so the "jet" vs. "corona" models are best in- 
terpreted as statements on which component is more dom- 
inant. Furthermore, factors other than just luminosity are 
important, as evidenced, e.g., by XTE J1550-564 showing 
different contributions of observed optically thin jet emis- 
sion to the X-ray (at similar luminosities) depending on if 
the black h ole is heading into ou tburst or declining back into 
quiescence (|Russell et al.ll201ol '). 

The above complications make using the FP to distin- 
guish X-ray processes a challenging, but not impossible, ex- 
ercise. That inclusion of accreting black holes at only certain 
accretion rates affects the intrinsic scatter and inferred co- 
efficients supports that GBH accretion states can indeed be 



mapped to different classes of AGN; thus, AGN unification 
is clearly more comp licated than just orientation (see, e.g., 
[Richards et al.|[201ll ). Looking toward the future, further ex- 
ploration quantifying how parameter space beyond accretion 
rate, for example, environment, black hole spin, etc., can af- 
fect GBH-AGN mappings would be very exciting. 
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